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ABSTRACT 

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, 
hard X-ray, and y-ray bands, and we compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency 
observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique 
data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost 
all the BL Lac objects have been detected by Fermi Large Area Telescope (LAT), whereas 30 to 40% of the flat-spectrum radio quasars (FSRQs) in 
the radio, soft X-ray, and hard X-ray selected samples are still below the y-ray detection limit even after integrating 27 months of Fermi-LAT data. 

The radio to sub-millimetre spectral slope of blazars is quite flat, with (a) ~ 0 up to about 70 GHz, above which it steepens to (a) 0.65. BL 

Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (v 8 cak ) in 
the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples with < v 8 cak ) = 10 13 1±0 1 Hz, while the mean inverse-Compton 
peak frequency, (yj 6: ak ), ranges from 10 21 to 10 22 Hz. The distributions of v 8 c . ik and of vj^ ak of BL Lacs are much broader and are shifted to higher 
energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars ranges from less than 
0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, 
with y-ray selected blazars peaking at ~ 7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption 
that the blazar power budget is largely dominated by high-energy emission is a selection effect. A comparison of our multi-frequency data with 
theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the y-ray detected blazars 
in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the 
correlation between bolometric luminosity and v 8 eak predicted by the blazar sequence. 


1. Introduction 


Blazars are jet-dominated extragalactic objects characterized by 
the emission of strongly variable and polarized non-thermal ra- 
diation across the entire electromagnetic spectrum, from radio 
waves to y-rays (e.g., |Urry & Padovani[|1995j l. As the extreme 
properties of these sources are due to relativistic amplification 
of radiation emitted along a jet pointing very close to the line of 
sight (e.g., |Blandford & Rees|1978[[Urry & Padovanill995[ ), they 
are rare compared to objects pointing their jets at random angles 
and to radio quiet QSOs where the emitted radiation is due to 
thermal or reflection mechanisms ultimately powered by accre- 
tion onto a supermassive black hole (e.g., Abdo et al. 2010a| . 
Despite that, blazars’ strong emission at all wavelengths makes 
them the dominant type of extragalactic sources in the radio, 
microwave, y-ray, and TeV bands where accretion and other 
thermal emission processes do not produce significant amounts 
of radiation ([Toffolatti et al.|fl998[ |Giommi & Colafrancesco! 


2004 [Hartman et al. |1999j Abdo et al. 201 Oat Costamante & 


Ghisellini 2002[ Colafrancesco & Giommi|2006[ Weekes|2 008 ). 
For these reasons blazars are hard to find at optical and X-ray fre 
quencies, while they dominate the microwave and y-ray sky at 
high G alactic latitud es. The advent of the Fermi (Atwood et al. 


2009 ) and Planck [[] ( jTauber et al. 2010 Planck Collaboration 


201 1 a[ > satellites, that are surveying these two observing win- 
dows, combined with the versatility of the Swift observatory 
( Gehrels et al.|2004a| , is giving us the unprecedented opportu- 
nity to collect multi-frequency data for very large samples of 
blazars and open the way to a potentially much deeper under- 


1 Planck (http://www.esa.int/Planck) is a project of the European 
Space Agency - ESA - with instruments provided by two scientific 
Consortia funded by ESA member states (in particular the lead coun- 
tries: France and Italy) with contributions from NASA (USA), and tele- 
scope reflectors provided in a collaboration between ESA and a scien- 
tific Consortium led and funded by Denmark. 
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standing of the physics and demographics of these still puzzling 
objects. 

Blazars can be categorized by their optical properties and the 
shape of their broad-band spectral energy distributions (SEDs). 
Blazar SEDs always show two broad bumps in the log v - 
log vF v space; the lower-energy one is usually attributed to syn- 
chrotron radiation while the more energetic one is attributed to 
inverse-Compton scattering. Blazars showing strong and broad 
optical emission lines are usually called flat-spectrum radio 
quasars (FSRQs), while objects with no broad emission lines 
(i.e., rest frame equivalent width, EW, < 5 A) are called BL 
Lac objects. Padovani & Giommi ( 1995] ) introduced the terms 
LSP and HSP to distinguish between BL Lacs with low and 
high values of the peak frequency of the synchrotron bump 

( V peak)-^ 


Abdo et al. 


(2010a) extended this definition to all types 
of blazars and defined the terms LSP, ISP, and HSP (correspond- 
ing to low, intermediate, and high synchrotron peaked blazars) 
for the cases where v s . < 10 14 Hz, 10 14 Hz < v s . < 10 15 Hz, 

peak peak 

and Vp eak > 10 15 Hz, respectively. In the rest of this paper we will 
use the LSP/ISP/HSP nomenclature. 

It is widely recognized that one of the most effective ways 
of studying the physics of blazars is through the use of multi- 
frequency data, ideally simultaneous. There are several exam- 


pies of studies following this approach (e.g.. 

Giommi et al.|1995| 

von Montigny et al. 1995 

Sambruna et al. 

1996) Fossati et al. 

1998 

Giommi et al. 2002 

Nieppola et al. 2006; Padovani et al. 

2006 

i, but in most cases t 

te samples are heterogeneous and the 


data are sparse and non-simultaneous. 

Simultaneous and well sampled SEDs require the organi- 
zation of complex multifrequency observation campaigns, in- 
volving the coordination of observations from several observa- 
tories. Such large efforts have been carried out only rarely, and 
almost exclusively on the occasion of large flaring events in 


a few bright and well known blazars, e.g., 3C 454.3 (Giommi 
|et al.| |2006[ Abdo et al. 2 009 a [ Vercellone et al.| 2009) , 


Mkn421, (|Donnarumma et al.||2009[ |Abdo et al.||201 1]), and 
PKS 2155-304 ( |Aharonian et al'|2009p . 

Significant progress has been made recently with the publi- 
cation of a compilation of quasi-simultaneous SEDs of a large 
sample of y-ray bright blazars (Abdo et al. 2010a| l. This is an 
important step forward compared to previous compilations, as 
the sample presented is statistically representative of the popu- 
lation of bright y-ray selected blazars, and the data were quasi- 
simultaneous, that is collected within three months of the Fermi- 
LAT observations. 

With Planck , Swift, and Fermi - LAT simultaneously in or- 
bit, complemented by other space and ground-based observato- 
ries, it is now possible to assemble high-quality multi-frequency 
datasets that allow us to build simultaneous and well sampled 
broad-band spectra of large and statistically well-defined sam- 
ples of AGNs. 

In this and a companion paper ( Planck Collaboration|201 le} 
we present the first results of a large cooperative program be- 
tween the Planck, Fermi- LAT, and Swift satellites and a num- 
ber of ground-based observatories, carried out to collect multi- 
frequency data on large samples of blazars selected using differ- 
ent criteria and observed when the sources lie in the field of view 
of Planck. 

In this paper we concentrate on blazars selected in the soft 
X-ray, hard X-ray, and y-ray bands. We present the simultaneous 
data, test for flux correlations, and estimate some key parame- 
ters characterizing the SEDs. We then compare the results ob- 
tained for the different samples. Detailed fits to models, variabil- 


ity studies, and more complete theoretical interpretations will be 
presented elsewhere. 

Throughout this paper we define the radio-to-submillimetre 
spectral index a by S (v) oc v", and we adopt a ACDM cosmol- 
ogy with //(> = 70kms~ 1 Mpc~', D„, = 0.27, and Da = 0.73 
(Komatsu et a l.|201l) l. 


2. Sample selection 

In order to explore the blazar parameter space from different 
viewpoints we have used several different criteria to select the 
blazars to be observed simultaneously by Planck, Swift, and 
Fermi. In this paper we consider three samples of blazars that 
are flux-limited in the high-energy part of the electromagnetic 
spectrum: soft X-ray (0. 1-2.4 keV) from the ROSAT All-Sky 
Survey Bright Source Catalog (1RXS, Voges et ak| |l999 here- 
after referred to as the RASS sample), hard X-ray (15-150 keV) 


from the Swift -BAT 54 month source catalog (Cusumano et al. 


2010 hereafter referred to as the BAT sample), and y-ray from 
the Fermi-b AT 3-month Bright AGN Source List (Abdo et al. 


2009b hereafter referred to as the Fermi-LAT sample). 

These high-energy-selected samples are complemented by 
a radio flux-limited sample of northern sources (hereafter re- 
ferred to as the radio sample) which is presented in the com- 
panion paper (Planck Collaboration| [2011e i. We use these four 
samples, defined in widely different parts of the electromagnetic 
spectrum, to try to disentangle the intrinsic properties of blazars 
from the heavy selection effects that often afflict blazar samples. 
Altogether we consider 175 sources. 

For the classification as different blazar types we rely on the 
Roma-BZCAT catalog (Massaro et al. 2010[ ) which is a com- 
pilation of known blazars that have been carefully checked in 
order to determine their blazar type in a uniform and reliable 
way. Massaro et al. ( 2010[ ) divide blazars into three categories: 
BZQ/FSRQ, in which the optical spectrum shows broad emis- 
sion lines; BZB/BL Lac objects, in which the optical spectrum is 
featureless or shows only absorption lines from the host galaxy; 
BZU/Uncertain type, comprising objects for which the authors 
could not find sufficient data to safely determine the source 
classification, and objects that show peculiar characteristics (see 
Massaro et al.|2010 for details). According to this classification 


96 of our objects are of the FSRQ type, 40 are BL Lacs, and 
the rest are of uncertain type. About 160 have been observed 
by Swift simultaneously with Planck, mostly through dedicated 
target of opportunity (ToO) pointings. In the following we de- 
scribe the selection criteria for each high-energy selected sam- 
ple. Details for the radio selected sample are given in Planck 
|Collaboration| ( |201 lej ). 

2. 1 . The issue of blazar classification 

This classification of blazars into featureless (BL Lacs) and 
broad-lined objects (FSRQs), although very simple in principle, 
is not unambiguous, nor it is robust. The borderline between the 
two blazar subclasses, namely 5 A in the source rest frame for 
the EW of any emission line, was originally defined as a result of 
the optical identification campaigns of the sources discovered in 
the first well defined and complete samples of (bright) radio and 
X-ray selected objects ( Stickel et al.||199l)|Stocke et al.|1991) . 
However, we now know that well-known BL Lac objects like 
OJ 287 - and BL Lac itself - exhibit emission lines with EW well 


above the 5 A limit on some occasions (Vermeulen et al. 1995 


Corbett et al. 1996) . Several other BL Lac objects show power- 
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ful emission lines with EW just below, and sometimes above the 
5 A threshold, depending on the variable continuum level (see 
e.g. Lawrence et a l ,| 1 996[ , [Ghi se Mini et al. 201 1[ and references 
therein). Well-known FSRQs like 3C 279 also appear nearly fea- 
tureless during bright states (Pian et al.~| |1999 1. The recent de- 
tection of broad Lyman-tr emission in the UV spectrum of clas- 
sical BL Lacs like Mkn421 and Mkn501 ( |Stocke et al.||201 T] > 
further contributes to the blurring of the distinction between the 
two types of blazars. 

The boundary between BL Lac objects and radio galaxies 
is also fragile as BL Lacs are distinguished from radio galax- 
ies when the 4000 A Ca H&K break (a stellar absorption fea- 
ture in the host galaxy) is diluted by non-thermal radiation more 
than a certain amount first quantified b y|Stocke et al.|(|1991| l and 
then revised by Marcha et al. ( 1996 1, Landt et al. ( |2002| i and 
Landt et al. ( 2004| > , The level of non-thermal blazar light around 
4000 A reflects the intrinsic radio power of the jet; it can be 
highly variable and depends strongly on the position of the peak 
of the synchrotron emission thus making the border between BL 
Lacs and radio galaxies quite uncertain. 

In this work we relied on the blazar classification given in the 
Roma-BZCAT catalog (Massaro et al. |2010 i which re-assessed 
the blazar subclass of each object after a critical review of the op- 
tical data available in the literature and in large public databases 
like the SDSS ( [York et al.|2000| . Despite that, some uncertain- 
ties remain, and our conclusions about the differences between 
BL Lacs and FSRQs could suffer from that. However, the large 
size of our samples ensures that a few misclassifications should 
not significantly affect our results. To assess the impact of blazar 
misclassification and of transitional objects in a quantitative way 
it is necessary to perform detailed simulations. 


2.2. The Fermi-LAT (y-ray flux limited) sample 

The y-ray flux limited sample has_been created from the 
Fermi LAT Bright Source Lis|^|( |Abdo et al.|2009b i. We have se- 


lected all the high Galactic latitude (\b\ > 10°) blazars detected 
with high significance (TS > 100/^] To further reduce the size 
of the sample and to ensure that all the sources are well above 
the Planck sensitivity limit for one survey, we have considered 
only the sources with radio flux density (taken from BZCAT) 
S 5 ghz > 1 Jy. We realize that this is a double cut, with a TS 
limit at y-ray energies and a flux-density limit in the radio band. 
Since a TS limit translates into different y-ray flux limits de- 
pending on the y-ray spectral slope, with higher sensitivity to 
flat-spectrum sources (see fig. 7 of Abdo et aT][2009b] l, for our 
statistical considerations we will also consider the subsample 
with a flux cut of F(E > 100 MeV) > 8 x 1(L 8 phcm _2 s _I , 
which removes such a dependence. The sample so defined in- 
cludes 50 sources, 40 of which are brighter than the y-ray flux 
limit of 8 x 1 0 8 ph ctrT 2 3 s 1 . Relaxing the radio flux density cut 
would give a purely y-ray flux limited sample and increase the 
number of sources to ~ 70, but with about 40-50% of the objects 
with S 5 ghz < 1 Jy undetected by Planck. We will consider this 
exclusively y-ray flux limited sample in the future when more 
Swift simultaneous observations will be available. 

Details are presented in Table [T] where column 1 gives the 
source common name, column 2 gives the Fermi-LAT name as 
it appears in Abdo et al. (2010b[), columns 3 and 4 give the 


2 http : / /www . asdc . asi . it/f ermibsl/ 

3 The Test Statistic (TS) is defined as TS = -21n(Lo/Li) with Lq 
the likelihood of the null-hypothesis model and the likelihood of a 
competitive model (see e.g. Abdo et al. 2010c). 


source position in equatorial coordinates, columns 5, 6, and 7 
give the redshift, V magnitude, and X-ray flux (0. 1-2.4 keV) 
from BZCAT ( |Massaro et al.|20K)| ), column 8 gives the 1 .4 GHz 
or 843 MHz flux density from NVSS ( jCondon et al,||1998) > or 
from SUMSS ( |Mauch et al.||2003|) with Dec < -40°), column 

9 gives the y-ray flux from Abdo et al. ( 2010b| ) p| and column 

10 gives the date of the Swift ToO observation made when the 
source was within the Planck field of view. 


2.3. The Swift/BAT (hard X-ray flux limited) sample 


We have defined our hard X-ray flux limited sample starting 
from the Swift-BAT 54 month source catalog (Cusumano et al. 
20 1 0 ), and selecting all the sources identified with blazars with 


X-ray flux > 11L 1 1 erg cnT 2 s 1 in the 15-150 keV energy band. 
The BAT catalog includes 70 known blazars that satisfy the X- 
ray flux condition, but many of them are too faint to be detected 
at millimetre wavelengths by Planck. Therefore, although a pure 
X-ray selection would be preferable, we have decided to add 
a mild radio flux-density constraint (Ssghz > 100 mJy, with 
S 5 ghz taken from BZCAT) to select only sources that can be 
detected by Planck or for which Planck will be able to provide 
meaningful upper limits, leaving enough sources to build a sta- 
tistically sizable sample. The list of the 34 sources included in 
this sample is given in Table [2] The column description is the 
same as for Table Q] 


2.4. The ROSAT/RASS (soft X-ray flux limited) sample 

The soft X-ray flux limited sample was defined starting from the 
RASS catalog (1RXS) (Voges et al. |1999| >, and selecting all the 
sources identified with blazars with count rate larger than 0.3 
counts/s in the 0. 1-2.4 keV energy band, and radio flux density 
from BZCAT S 5 ghz > 200 mJy. The reasons to use a mild radio 
flux constraint are the same as for the hard X-ray flux limited 
sample; for this case, however, we have chosen 200 mJy to fur- 
ther reduce the size of the sample to be comparable with the y- 
ray and hard X-ray samples. The list of the 43 sources included 
in this sample is given in Table [3] The column description is the 
same as for Tables |T] and [2] 


2.5. The radio flux-density limited sample 


The radio flux-density limited sample is presented in the com- 
panion Planck paper (Planck Collaboration 201 le| > where all 
the observational detai s are given. The sample consists of 104 
bright northern and equatorial radio-loud AGN characterized by 
S 37 ghz > 1 Jy as measured with the Metsahovi radio telescope. 


Although the samples are defined by different criteria, four 
sources are common to all samples. These are the well-known 
objects 3C273, 3C279, Mkn501 and 3C 454.3, which are 
among the brightest objects across the entire electromagnetic 
spectrum. A summary of the number of sources common to more 
than one sample is given in Table [4] 


4 We give average y-ray fluxes from the 1-year Fermi catalog rather 
than the three-month fluxes that were used to define the sample for con- 
sistency with Table s|2] and [3] 

5 http://v7ww.asdc.asi.it/bat54/ 
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Table 1. The Fermi-LAT (y-ray TS/flux limited) sample. 


Source name 

Fermi-LAT name (1FGL) 

R.A. 

(J2000) 

Dec. 

(J2000) 

z 

Rmag 

X-ray Flux 
0. 1-2.4 keV 

b 

Flux Density 
1.4 GHz 
mJy 

Fermi Flux 
1-100 GeV 

Swift obs. date 

Blazar type 

lJy 01 18-272 

1FGLJ 0120.5-2700 

01 20 31.6 

-27 01 24 

0.557 

15.5 

0.72 

934 

3.7 ± 0.4 

d 

BL Lac - LSP 

S40133+47 

1FGLJ 0137.0+4751 

01 36 58.5 

47 51 29 

0.859 

17.6 

1.04 

1138 

9.6 ± 0.6 

2010-02-05 

FSRQ -LSP 

PKS 0202-17 

1FGLJ 0205.0- 1702 

02 04 57.6 

-17 01 18 

1.740 

17.3 

0.57 

1220 

1.5 ±0.3 

2010-01-08 

FSRQ - LSP 

PKS 0208-512 

1FGLJ 0210.6-5 101 

02 10 46.2 

-51 01 01 

1.003 

14.8 

0.75 

3493 

7.1 ±0.6 

2009-11-26 

Uncertain - LSP 

PKS 0215+015 

1FGLJ 0217.9+0144 

02 17 48.9 

01 44 49 

1.715 

18.7 

2.56 

751 

6.0 ± 0.5 

d 

FSRQ - LSP 

lJy 0218+357 

1FGLJ 022 1.0+3555 

02 21 05.5 

35 56 14 

0.944 

20.0 

0.85 

1707 

6.4 ± 0.5 

2010-08-19 

Uncertain - LSP 

4C 28.07 

1FGLJ 0237.9+2848 

02 37 52.4 

28 48 09 

1.213 

18.8 

0.58 

2197 

3.7 ± 0.4 

2010-02-05 

FSRQ - LSP 

PKS 0235+164 

1FGLJ 0238.6+1637 

02 38 38.9 

16 36 59 

0.940 

18.5 

1.24 

1941 

32.7 ± 1.1 

2010-01-30 

BL Lac - LSP 

PKS 0332-403 

1FGLJ 0334.2-4010 

03 34 13.6 

-40 08 25 


17.5 

0.73 

1042 

3.8 ± 0.4 

2010-01-17 

BLLac -LSP 

PKS 0420-01 

1FGLJ 0423.2-01 18 

04 23 15.8 

-01 20 33 

0.916 

16.7 

1.39 

2726 

5.6 ± 0.5 

2009-08-27 

FSRQ - LSP 

PKS 0426-380 

1FGLJ 0428.6-3756 

04 28 40.4 

-37 56 19 

1.030 

16.3 

0.42 

753 

25.7 ± 1.0 

2010-08-17 

BL Lac - LSP 

PKS 0454-234 

1FGLJ 0457.0-2325 

04 57 03.1 

-23 24 52 

1.003 

17.9 

0.91 

1727 

32.5 ± 1.1 

2010-02-25 

FSRQ - LSP 

PKS 0528+134 

1FGLJ 053 1.0+ 1331 

05 30 56.4 

13 31 55 

2.070 

18.9 

0.80 

1556 

4.0 ± 0.5 

2009-09-24 

FSRQ - LSP 

PKS 0537-441 

1FGLJ 0538.8-4404 

05 38 50.3 

-44 05 08 

0.892 

16.0 

2.10 

3729 

21.3 ±0.9 

2010-03-03 

BLLac -LSP 

PKS 0735+17 

1FGLJ 0738.2+1741 

07 38 07.3 

17 42 19 

0.424 

14.5 

0.97 

2258 

4.4 ± 0.5 

2010-10-07 

BLLac -ISP 

S4 08 14+425 

1FGLJ 08 18.2+4222 

08 18 16.0 

42 22 45 

0.530 

19.6 

0.32 

1091 

8.7 ± 0.6 

2010-10-15 

BL Lac - LSP 

OJ 535 

1FGLJ 0825.0+5555 

08 24 47.2 

55 52 42 

1.417 

18.1 

0.64 

1449 

0.9 ± 0.3 

2010-03-28 

FSRQ - LSP 

PKS 0851+202 

1FGLJ 0854.8+2006 

08 54 48.8 

20 06 30 

0.306 

14.4 

1.72 

1512 

2.7 ± 0.4 

2010-04-10 

BL Lac - LSP 

S4 09 17 +44 

1FGLJ 0920.9+4441 

09 20 58.4 

44 41 54 

2.190 

19.2 

1.04 

1017 

14.0 ± 0.7 

2009-10-29 

FSRQ - LSP 

4C55.17 

1FGLJ 0957.7+5523 

09 57 38.1 

55 22 57 

0.896 

16.8 

0.51 

3079 

10.5 ± 0.6 

2009-11-01 

FSRQ - LSP 

4C 01.28 

1FGLJ 1058.4+0134 

10 58 29.6 

01 33 58 

0.888 

17.6 

1.08 

3220 

7.1 ±0.6 

2009-12-03 

Uncertain - LSP 

PKS 1057-79 

1FGLJ 1058.1-8006 

10 58 43.3 

-80 03 54 

0.581 

17.3 

0.43 

534 

2.2 ± 0.4 

2010-08-30 

BLLac -LSP 

PKS 1127-145 

1FGLJ 1130.2-1447 

11 30 07.0 

-14 49 27 

1.184 

16.0 

1.39 

5622 

2.4 ± 0.4 

2009-12-28 

FSRQ - LSP 

PKS 1144-379 

1FGLJ 1146.9-3812 

11 47 01.3 

-38 12 11 

1.048 

15.7 

0.91 

1804 

2.4 ± 0.4 

2010-06-24 

BLLac -LSP 

4C 29.45 

1FGLJ 1159.4+2914 

11 5931.8 

29 14 44 

0.729 

16.4 

0.84 

2031 

5.3 ± 0.5 

2010-05-28 

FSRQ - LSP 

ON 231 

1FGLJ 1221.5+2814 

1221 31.6 

28 13 58 

0.102 

14.3 

1.30 

732 

6.9 ± 0.5 

2009-12-10 

BLLac -ISP 

3C273 

1FGLJ 1229.1+0203 

12 29 06.7 

02 03 08 

0.158 

14.1 

63.11 

54991 

9.6 ± 0.6 

2010-01-11 

FSRQ - LSP 

PKS 1244-255 

1FGLJ 1246.7-2545 

12 46 46.8 

-25 47 49 

0.635 

16.7 

1.30 

1165 

8.1 ±0.6 

2010-01-25 

FSRQ - LSP 

3C279 

1FGLJ 1256.2-0547 

1256 11.1 

-05 47 21 

0.536 

15.0 

20.90 

9711 

32.4 ± 1.1 

2010-01-15 

FSRQ - LSP 

lJy 1308+326 

1FGLJ 1310.6+3222 

13 10 28.6 

32 20 43 

0.997 

19.6 

0.53 

1687 

6.8 ± 0.5 

2009-12-12 

Uncertain - LSP 

PKS 1502+106 

1FGLJ 1504.4+1029 

15 04 24.9 

10 29 39 

1.839 

19.5 

0.16 

1774 

67.0 ± 1.6 

2010-07-29 

FSRQ - LSP 

4C -05.64 

1FGLJ 151 1.1-0545 

15 10 53.5 

-05 43 07 

1.191 

16.9 

1.16 

3569 

2.1 ±0.4 

d 

FSRQ - LSP 

APLib 

1FGLJ 1517.8-2423 

15 17 41.8 

-24 22 19 

0.048 

12.6 

1.05 

2042 

5.6 ± 0.5 

2010-02-20 

BL Lac - LSP 

OS -237.8 

1FGLJ 1625.7-2524 

16 25 46.8 

-25 27 38 

0.786 

20.6 


2521 

4.4 ± 0.6 

2010-08-26 

Uncertain - LSP 

4C 38.41 

1FGLJ 1635.0+3808 

16 35 15.4 

38 08 04 

1.814 

17.3 

0.17 

2726 

6.8 ± 0.5 

2010-03-07 

FSRQ - LSP 

NRA0512 

1FGLJ 1642.5+3947 

16 40 29.6 

39 46 46 

1.660 

17.5 

0.15 

976 

5.6 ± 0.5 

2010-08-06 

FSRQ - LSP 

Mkn 501 

1FGLJ 1653.9+3945 

16 53 52.2 

39 45 36 

0.033 

8.3 

36.90 

1558 

8.3 ±0.6 

2010-03-21 

BLLac-HSP 

OT081 

1FGLJ 1751.5+0937 

17 51 32.8 

09 39 00 

0.322 

17.0 

1.18 

623 

6.4 ± 0.6 

2010-04-01 

Uncertain - LSP 

S5 1803+784 

1FGLJ 1800.4+7827 

18 00 45.6 

78 28 04 

0.680 

14.7 

0.79 

2223 

3.0 ± 0.4 

2009-10-13 

BLLac -LSP 

2E 1908.2-2011 

1FGLJ 1911.2-2007 

19 11 09.6 

-20 06 55 

1.119 

18.9 

1.77 

2714 

4.5 ± 0.5 

2009-10-04 

FSRQ - LSP 

PMNJ 1923-2104 

1FGLJ 1923.5-2104 

19 23 32.1 

-21 04 33 

0.874 

16.6 

0.77 

3167 

11.9 ±0.7 

2010-09-30 

FSRQ - LSP 

lJy 2005-489 

1FGLJ 2009.5-4849 

20 09 25.3 

-48 49 53 

0.071 

11.0 

33.24 

1282 

5.0 ± 0.5 

2009-10-05 

BLLac-HSP 

PKS 2052-47 

1FGLJ 2056.3-47 14 

20 56 16.3 

-47 14 47 

1.491 

18.3 

0.56 

2223 

4.6 ± 0.5 

2010-10-18 

FSRQ - LSP 

S3 2141+17 

1FGLJ 2143.4+1742 

21 43 35.5 

17 43 48 

0.213 

14.4 

0.63 

651 

4.9 ± 0.5 

2009-11-20 

FSRQ - LSP 

lJy 2144+092 

1FGLJ 2147.2+0929 

21 47 10.1 

09 29 46 

1.113 

16.9 

1.49 

934 

4.1 ±0.4 

d 

FSRQ - LSP 

BLLac 

1FGLJ 2202.8+4216 

22 02 43.2 

42 16 40 

0.069 

12.5 

1.58 

6051 

7.1 ±0.6 

2009-12-23 

BLLac -LSP 

PKS 2204-54 

1FGLJ 2207.8-5344 

22 07 43.7 

-53 46 33 

1.215 

18.2 

0.52 

1526 

1.6 ±0.3 

2010-05-03 

FSRQ - LSP 

PKS 2227-08 

1FGLJ 2229.7-0832 

22 29 40.0 

-08 32 54 

1.560 

16.8 

8.74 

968 

4.6 ± 0.5 

2009-11-19 

FSRQ - LSP 

4C 11.69 

1FGLJ 2232.5+1 144 

22 32 36.4 

1 1 43 50 

1.037 

16.5 

1.26 

7202 

4.1 ±0.4 

2009-11-29 

FSRQ - LSP 

3C 454.3 

1FGLJ 2253.9+1608 

22 53 57.7 

16 08 53 

0.859 

13.2 

7.80 

12657 

46.2 ± 1.3 

2009-12-14 

FSRQ - LSP 


Notes. 843 MHz flux is reported for sources with Dec < -40°; (6) Units of 10 12 ergcm 2 s *; (c) Units of 10 10 phcm 2 s 1 
(rf) Swift simultaneous observation not available. w Optical spectrum completely featureless or not available, redshift unknown 


3. Data analysis 

3. 1 . Ground-based follow-up data 


version 30-0ct-2009. Uranus was used as a calibrator 
for the flux densities. 


Following the launch of Planck, several follow-up programs 
with ground-based facilities started collecting simultaneous ra- 
dio and optical data. In this paper we use data from the observa- 
tories listed in Table 0 


3.1.1. APEX 

Some sources from our sample were observed in the submil- 
limetre domain with the 12-m Atacama Pathfinder Experiment 
(APEX) in Chile. The observations were made using the 
LABOCA bolometer array centered at 345 GHz. Data were 
taken at two epochs: September 3-4, 2009, and November 12, 
2009. The data were reduced using the script package mini- 


3.1.2. ATCA-PACO 


The Planck-AICA Co-eval Observations (PACO) project 
dMassardi et al.|[201 1 a| |Bonavera et al.||2011| ) observed 480 
sources selected from the Australia Telescope 20 GHz catalogue 
(AT20G, Massardi et al. 201 lb| ), with the Australia Telescope 
Compact Array (ATCA) in the frequency range 4.5-40 GHz, at 
several epochs close in time to the Planck observations in the pe- 
riod July 2009 to August 2010. The PACO sample is a complete, 
flux-density limited, and spectrally selected sample of southern 
sources, with the exception of the region with Galactic latitude 


http : //www . submm . caltech . edu/ - share/ crush/ 
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Table 2. The Swift-BAI (hard X-ray flux limited) sample. 


Source name 

Fermi - LAT name 

R.A. 

(J2000) 

Dec. 

(J2000) 

z 

Rmag 

X-ray Flux 
0. 1-2.4 keV 

b 

Flux Density 
1.4 GHz 
mJy 

Fermi-L AT Flux 
1-100 GeV 

Swift obs. date 

Blazar type 

III ZW 2 


00 10 31.0 

10 58 29 

0.089 

13.9 

6.14 

98 


2010-07-08 

FSRQ - LSP 

S5 0014+813 


00 17 08.4 

81 35 08 

3.387 

15.9 

0.77 

693 


2010-09-21 

FSRQ - ISP 

1ES 0033+595 

1FGLJ 0035.9+5951 

00 35 52.6 

59 50 03 

e 

17.2 

5.41 

147 

3.2 ± 0.5 

d 

BLLac -HSP 

Mkn 348 


00 48 47.1 

31 57 25 

0.015 

9.3 

0.12 

292 


d 

BLLac -LSP 

lJy 0212+735 

1FGLJ 0217.8+7353 

02 17 30.8 

73 49 32 

2.367 

18.8 

0.54 

2272 

1.0 ±0.4 

2010-09-11 

Uncertain - LSP 

NGC 1275 

1FGLJ 03 19.7+4130 

03 19 48.1 

41 30 42 

0.018 

12.3 

197.92 

22830 

17.3 ± 0.8 

2010-08-09 

FSRQ -LSP 

NR AO 140 

1FGLJ 0334.2+3233 

03 36 30.1 

32 18 29 

1.259 

16.6 

2.80 

2677 

1.0 ±0.4 

2010-08-25 

Uncertain - LSP 

3C120 


04 33 11.0 

05 21 15 

0.033 

13.8 

22.68 

3440 


2010-02-25 

Uncertain - LSP 

PKS 0521-36 

1FGLJ 0522.8-3632 

05 22 57.9 

-36 27 30 

0.055 

11.6 

10.42 

15620 

2.9 ± 0.4 

2010-03-05 

Uncertain - LSP 

PKS 0528+134 

1FGLJ 053 1.0+ 1331 

05 30 56.4 

13 31 55 

2.070 

18.9 

0.80 

1556 

4.0 ± 0.5 

2009-09-24 

FSRQ - LSP 

lJy 0537-286 

1FGLJ 0539. 1-2847 

05 39 54.2 

-28 39 55 

3.104 

19.0 

0.83 

862 

0.9 ± 0.0 

2010-03-12 

FSRQ - LSP 

PKS 0548-322 


05 50 40.6 

-32 16 20 

0.069 

13.1 

26.34 

344 


2010-03-12 

BLLac -HSP 

B2.2 0743+25 

1FGLJ 0746.6+2548 

07 46 25.8 

25 49 02 

2.979 

19.2 

0.38 

417 

0.7 ± 0.2 

2010-10-15 

FSRQ - LSP 

4C 71.07 

1FGLJ 0842.2+7054 

08 41 24.3 

70 53 42 

2.218 

16.8 

5.52 

3823 

1.2 ±0.3 

2010-03-21 

FSRQ - LSP 

Mkn 421 

1FGLJ 1104.4+3812 

11 04 27.3 

38 12 31 

0.030 

8.3 

180.94 

767 

26.1 ± 1.0 

2009-11-17 

BLLac -HSP 

PKS 1127-145 

1FGLJ 1130.2-1447 

11 30 07.0 

-14 49 27 

1.184 

16.0 

1.39 

5622 

2.4 ± 0.4 

2009-12-28 

FSRQ - LSP 

PKS 1219+04 

1FGLJ 1222.5+0415 

12 22 22.5 

04 13 15 

0.967 

17.1 

0.91 

800 

0.9 ± 0.3 

2010-07-17 

FSRQ - LSP 

3C273 

1FGLJ 1229.1+0203 

12 29 06.7 

02 03 08 

0.158 

14.1 

63.11 

54991 

9.6 ± 0.6 

2010-01-11 

FSRQ - LSP 

3C279 

1FGLJ 1256.2-0547 

12 56 11.1 

-05 47 21 

0.536 

15.0 

20.85 

9711 

32.4 ± 1.1 

2010-01-15 

FSRQ - LSP 

APLib 

1FGLJ 1517.8-2423 

15 1741.8 

-24 22 19 

0.048 

12.6 

1.05 

2042 

5.7 ± 0.5 

2010-02-20 

BLLac -LSP 

Mkn 501 

1FGLJ 1653.9+3945 

16 53 52.2 

39 45 36 

0.033 

8.3 

36.93 

1558 

8.3 ± 0.6 

2010-03-21 

BLLac -HSP 

ARP 102B 


17 19 14.4 

48 58 49 

0.024 

9.4 

0.59 

145 


2010-03-31 

Uncertain - ISP 

PKSB 1830-210 

1FGLJ 1833.6-2103 

18 33 39.8 

-21 03 39 

2.507 

16.6 

0.69 

10896 

10.7 ± 0.8 

2010-09-23 

FSRQ - LSP 

OV -236 

1FGLJ 1925.2-2919 

19 24 51.0 

-29 14 30 

0.352 

17.3 

2.42 

13387 

1.4 ±0.4 

2010-09-30 

FSRQ - LSP 

1ES 1959+650 

1FGLJ 2000.0+6508 

19 59 59.8 

65 08 54 

0.047 

11.9 

35.28 

250 

6.0 ± 0.5 

2009-09-26 

BLLac -HSP 

1 Jy 2126—158 


21 29 12.1 

-15 3841 

3.268 

16.5 

1.54 

590 


2010-05-03 

FSRQ - ISP 

4C 06.69 

1FGLJ 2148.5+0654 

21 48 05.4 

06 57 38 

0.999 

15.1 

1.46 

2589 

0.7 ± 0.3 

2009-11-21 

FSRQ - LSP 

PKS 2149-307 


21 51 55.5 

-30 27 53 

2.345 

17.4 

4.80 

1243 


2010-05-11 

FSRQ - LSP 

BLLac 

1FGLJ 2202.8+4216 

22 02 43.2 

42 16 40 

0.069 

12.5 

1.57 

6051 

7.1 ±0.6 

2009-12-23 

BLLac -LSP 

4C 31.63 


22 03 14.9 

31 45 38 

0.295 

14.3 

3.22 

2878 


2009-11-27 

FSRQ - LSP 

NGC 7213 


22 09 16.2 

-47 10 00 

0.006 

10.3 

35.34 

98 


2010-10-23 

Uncertain 

4C 11.69 

1FGLJ 2232.5+1 144 

22 32 36.4 

1 1 43 50 

1.037 

16.5 

1.26 

7202 

4.1 ±0.4 

2009-11-29 

FSRQ - LSP 

3C 454.3 

1FGLJ 2253.9+1608 

22 53 57.7 

16 08 53 

0.859 

13.2 

7.80 

12657 

46.2 ± 1.3 

2009-12-14 

FSRQ - LSP 

PKS 2325+093 

1FGLJ 2327.7+0943 

23 27 33.5 

09 40 09 

1.843 

18.8 

0.73 

741 

3.0 ± 0.4 

2010-07-18 

FSRQ - LSP 


Notes. (a> 843 MHz flux is reported for sources with Dec < -40°. (b) Units of 10 12 ergcm 2 s *; (c> Units of 10 10 phcm 2 s 1 
(d> Swift simultaneous observation not available. w Optical spectrum completely featureless or not available, redshift unknown 


\b\ < 5°. A total of 147 PACO point-like sources have at least 
one observation within 10 days of the Planck observations. 


3.1.3. Effelsberg and IRAM 

Quasi-simultaneous cm/mm radio spectra for a larger number 
of Planck blazars have been obtained within the framework 
of a Ferwimonitoring program of y-ray blazars (F-GAMMA: 
|Fuhrmann et al.| 2007 [ |Angelakis et al.|[2008| l on the Effelsberg 
100-m and IRAM 30-m telescopes. The frequency range was 
2.64-142 GHz. 

The Effelsberg measurements were conducted with the sec- 
ondary focus heterodyne receivers at 2.64, 4.85, 8.35, 10.45, 
14.60, 23.05, 32.00, and 43.00 GHz. The observations were 
performed quasi-simultaneously with cross-scans, that is slew- 
ing over the source position in azimuth and elevation with the 
number of sub-scans chosen to reach the desired sensitivity 
(for details, see Fuhrmann et al.||2008| Angelakis et al.||2008) . 
Pointing offset, gain, atmospheric opacity, and sensitivity cor- 
rections have been applied to the data. 

The IRAM 30-m observations were carried out with cali- 
brated cross-scans using the EMIR horizontal and vertical polar- 
ization receivers operating at 86.2 and 142.3 GHz. The opacity- 
corrected intensities were converted into the standard tempera- 
ture scale and finally corrected for small remaining pointing off- 
sets and systematic gain-elevation effects. Conversion to stan- 
dard flux density scale was based on frequent observations of 
primary calibrators (Mars, Uranus) and secondary calibrators 
(W3(OH), K3-50A, NGC 7027). 


From this program, radio spectra measured quasi- 
simultaneously with the Planck observations have been 
collected for a total of 37 Planck blazars during the period 
August 2009 to June 2010. 


3.1.4. Medicina 

Medicina Planck Experiment (SiMPlE, 
) is using the 32-m Medicina single dish to 
make almost simultaneous observations at 5 and 8.3 GHz of the 
263 sources of the NEWPS sample ( |Massardi et al.|2009) > with 
Dec > 0°, and partially overlapping with the PACO observations 
for -10° < Dec < 0°. The project began in June 2010 and fin- 
ished in June 2011, observing the sample several times through- 
out two complete Planck surveys. It does not overlap with the 
Planck first survey. 


The Simultaneous 
Procopio et al. |201 


3.1.5. Metsahovi 

The 37 GHz observations were made with the 13.7-m Metsahovi 
radio telescope using a 1 GHz bandwidth, dual-beam receiver 
centered at 36.8 GHz. The observations are ON-ON observa- 
tions, alternating the source and the sky in each feed horn. A 
typical integration time to obtain one flux density data point is 
1200-1400 s. The telescope detection limit at 37 GHz is ~0.2 Jy 
under optimal conditions. Data points with a signal-to-noise ra- 
tio below four are handled as non-detections. The flux-density 
scale is set by observations of DR 21. Sources NGC 7027, 
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Table 3. The ROSAT/RASS (soft X-ray flux limited) sample. 


Source name 

Fermi-hAT name 

R.A. 

(J2000) 

Dec. 

(J2000) 

z 

Rmag 

X-ray Flux 
0. 1-2.4 keV 

b 

Flux Density 
1.4 GHz 
mJy 

Fermi - LAT Flux 
1-100 GeV 

Swift obs. date 

Blazar type 

III ZW 2 


00 10 31.0 

10 58 29 

0.089 

13.9 

6.14 

98 


2010-07-08 

FSRQ - LSP 

GB6J 0214+5 145 


02 14 17.9 

51 44 52 

0.049 

16.5 

4.58 

294 


d 

BLLac-HSP 

3C120 


04 33 11.0 

05 21 15 

0.033 

13.8 

22.68 

3440 


2010-02-25 

Uncertain - LSP 

PKS 0521-36 

1FGLJ 0522.8-3632 

05 22 57.9 

-36 27 30 

0.055 

11.6 

10.42 

15620 

2.9 ± 0.4 

2010-03-05 

Uncertain - LSP 

PKS 0548-322 


05 50 40.6 

-32 16 20 

0.069 

13.1 

26.34 

344 


2010-03-12 

BLLac-HSP 

IRAS-L 06229-643 


06 23 07.6 

-64 36 20 

0.129 

13.7 

5.34 

274 


2010-08-18 

FSRQ - LSP 

4C 71.07 

1FGLJ 0842.2+7054 

08 41 24.3 

70 53 42 

2.218 

16.8 

5.52 

3823 

1.2 ±0.3 

2010-03-21 

FSRQ - LSP 

B2 0912+29 

1FGLJ 09 15.7+2931 

09 15 52.4 

29 33 24 

e 

15.0 

6.25 

342 

2.1 ±0.3 

2010-10-28 

BLLac -HSP 

PKS 0921-213 


09 23 38.8 

-21 35 47 

0.053 

12.8 

4.80 

268 


2010-05-02 

Uncertain - LSP 

1H 1013+498 

1FGLJ 1015.1+4927 

10 15 04.1 

49 26 00 

0.212 

15.1 

13.23 

378 

8.7 ± 0.6 

2010-04-24 

BLLac-HSP 

1RXSJ 105837.5+562816 

1FGLJ 1058.6+5628 

10 58 37.7 

56 28 1 1 

0.143 

14.0 

3.13 

228 

5.7 ± 0.5 

2010-04-18 

BLLac-HSP 

PKS 1124-186 

1FGLJ 1126.8-1854 

11 27 04.3 

-18 57 17 

1.048 

19.2 

5.33 

536 

2.4 ± 0.4 

2010-06-10 

FSRQ - LSP 

B2 1128+31 


11 31 09.4 

31 14 05 

0.289 

15.8 

5.02 

370 


2009-11-28 

FSRQ - HSP 

S5 1133+704 

1FGLJ 1136.6+7009 

1 1 36 26.4 

70 09 27 

0.045 

11.0 

35.08 

327 

1.7 ± 0.3 

2009-10-27 

BLLac-HSP 

4C 49.22 


11 53 24.4 

49 31 08 

0.334 

16.9 

3.31 

1572 


200911-17 

FSRQ - LSP 

ON 325 

1FGLJ 1217.7+3007 

12 17 52.0 

30 07 00 

0.130 

14.5 

24.90 

572 

6.7 ± 0.6 

2009-12-03 

BLLac-HSP 

PKS 1217+02 


12 20 11.8 

02 03 42 

0.241 

15.6 

2.78 

672 


2010-06-24 

FSRQ - LSP 

3C273 

1FGLJ 1229.1+0203 

12 29 06.7 

02 03 08 

0.158 

14.1 

63.11 

54991 

9.6 ± 0.6 

2010-01-11 

FSRQ - LSP 

PG 1246+586 

1FGLJ 1248.2+5820 

12 48 18.7 

58 20 28 

e 

14.5 

3.99 

245 

4.5 ± 0.4 

2010-05-20 

BLLac -ISP 

3C279 

1FGLJ 1256.2-0547 

1256 11.1 

-05 47 21 

0.536 

15.0 

20.85 

9711 

32.4 ± 1.1 

2010-01-15 

FSRQ - LSP 

lJy 1302-102 


13 05 33.0 

-1033 19 

0.286 

14.4 

4.20 

711 


d 

FSRQ - ISP 

GB6B 1347+0955 


13 50 22.1 

09 40 10 

0.133 

13.6 

3.74 

293 


2010-07-18 

Uncertain - ISP 

1WGAJ 1407.5-2700 


14 07 29.7 

-27 01 04 

0.022 

9.7 

15.28 

646 


2010-02-12 

Uncertain - HSP 

3C 298.0 


14 19 08.1 

06 28 34 

1.437 

16.4 

0.00 

6100 


2010-07-30 

Radio Galaxy 

BZQJ 1423+5055 


14 23 14.1 

50 55 37 

0.286 

15.1 

3.35 

178 


2010-07-13 

FSRQ - HSP 

PG 1424+240 

1FGLJ 1426.9+2347 

14 27 00.3 

23 48 00 

e 

14.5 

3.57 

430 

10.2 ± 0.6 

2010-01-22 

BLLac -ISP 

1RXSJ 145603.4+504825 


14 56 03.6 

50 48 25 

0.479 

18.1 

13.02 

220 


2009-12-25 

BLLac-HSP 

BZQJ 1507+0415 


15 07 59.7 

04 15 12 

1.701 

19.0 

6.11 

167 


2010-08-05 

FSRQ - LSP 

PG 1553+1 13 


15 55 43.0 

1 1 1 1 24 

e 

14.0 

17.85 

312 


2010-02-05 

BLLac -ISP 

WE 1601+16W3 


16 03 38.0 

15 54 02 

0.110 

13.4 

4.14 

97 


2010-08-14 

Uncertain - HSP 

3C345 

1FGLJ 1642.5+3947 

16 42 58.8 

39 48 37 

0.593 

15.0 

2.52 

7099 

5.6 ± 0.5 

2010-03-06 

FSRQ - LSP 

Mkn 501 

1FGLJ 1653.9+3945 

16 53 52.2 

39 45 36 

0.033 

8.3 

36.93 

1558 

8.3 ± 0.6 

2010-03-21 

BLLac-HSP 

1ES 1741 + 196 

1FGLJ 1744.2+1934 

17 43 57.8 

19 35 09 

0.084 

12.7 

4.23 

301 

1.1 ±0.3 

d 

BLLac -ISP 

PKS 1833-77 


18 40 38.4 

-77 09 28 

0.018 

8.3 

5.93 

1108 


2010-03-11 

Uncertain - ISP 

1ES 1959+650 

1FGLJ 2000.0+6508 

19 59 59.8 

65 08 54 

0.047 

11.9 

35.28 

250 

6.0 ± 0.5 

2009-09-26 

BLLac-HSP 

lJy 2005-489 

1FGLJ 2009.5-4849 

20 09 25.3 

-48 49 53 

0.071 

11.0 

33.24 

1282 

5.0 ± 0.5 

2009-10-05 

BLLac-HSP 

PKS 2149-307 


21 51 55.5 

-30 27 53 

2.345 

17.4 

4.80 

1243 


2010-05-11 

FSRQ - LSP 

NGC7213 


22 09 16.2 

-47 10 00 

0.006 

10.3 

35.34 

98 


2010-10-23 

Uncertain - HSP 

PKS 2227-399 


22 30 40.2 

-39 42 52 

0.318 

16.0 

4.23 

369 


2010-05-09 

Uncertain - ISP 

3C 454.3 

1FGLJ 2253.9+1608 

22 53 57.7 

16 08 53 

0.859 

13.2 

7.80 

12657 

46.2 ± 1.3 

2009-12-14 

FSRQ - LSP 

PKS 2300-18 


23 03 02.9 

-18 41 25 

0.129 

15.5 

5.16 

861 


2010-05-30 

Uncertain - ISP 

PKS 2331-240 


23 33 55.2 

-23 43 40 

0.048 

11.4 

3.92 

782 


2010-06-05 

Uncertain - ISP 

1ES 2344+514 

1FGLJ 2347. 1+5 142 

23 47 04.8 

51 42 17 

0.044 

10.7 

7.71 

250 

1.4 ±0.3 

2010-01-17 

BLLac-HSP 


Notes. (a) 843 MHz flux is reported for sources with Dec < -40°. (i) Units of 10 12 erg cm 2 s 1 ; (c) Units of 10 10 ph cm 2 s 1 
(rf) Swift simultaneous observation not available. Ic) Optical spectrum completely featureless or not available, redshift unknown 


Table 4. Summary of the samples, blazar types and selection methods considered in this paper. 


Sample 

Selection 

band 

No. of 
sources 

Blazars 

FS/BL/Unc. 

Other AGN 

RASS 

Sources in common 
BAT Fermi - LAT Radio 

All 

RASS 

Soft X-ray 

43 

15/16/11 

1 


12 

5 

9 

4 

BAT 

Hard X-ray 

34 

21/7/6 


12 


9 

16 

4 

Fermi - LAT * 

y-ray 

50 

28/16/6 


5 

9 


23 

4 

Fermi-LAT FL” 

y- ray 

40 

27/8/5 


3 

7 


19 

3 

Total this paper 


105 

52/32/20 

1 






Radio 

radio 

104 

73/18/10 

3 

9 

16 

23 


4 


* Total Fermi - LAT sample (TS limited), 

** Flux limited Fermi - LAT sample F(E > 100 Me V) > 8 x 10 -8 ph cnT 2 s _1 


3C274, and 3C84 are used as secondary calibrators. A de- 
tailed description of the data reduction and analysis is given in 
(Terasranta et aT7| |1998] >. The error estimate in the flux density 
includes the contribution from the measurement rms and the un- 
certainty of the absolute calibration. 


3.1.6. OVRO 


Some of the sources in our samples were monitored at 15 GHz 
using the 40-metre Telescope of the Owens Valley Radio 
Observatory as part of a larger monitoring program ( | Richards] 
et al. 12011 1. The flux density of each source was measured ap- 


proximately twice weekly, with occasional gaps due to poor 
weather or instrumental problems. The telescope is equipped 
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Table 5. Optical and radio observatories participating in the Planck multi-frequency campaigns. 


Radio observatory 

Frequencies [GHz] 

APEX, Chile 

345 


4.7, 5.2, 5.8, 6.3, 8.2, 8.7, 9.3, 9.8, 17.2, 17.7, 18.3, 18.8 


23.2, 23.7, 24.3, 24.8, 32.2, 32.7, 33.3, 33.8, 38.2, 38.7, 39.3, 39.8 

Effelsberg, Germany 

2.64, 4.85, 8.35, 10.45, 14.60, 23.05, 32.00, 43.00 

IRAM, Spain 

86.2, 142.3 

Medicina, Italy 

5, 8 

Metsahovi, Finland 

37 

OVRO, USA 

15 

RATAN, Russia 

1.1, 2.3, 4.8, 7.7, 11.2,21.7 

UMRAO, USA 

4.8, 8.0, 14.5 

VLA, USA 

5, 8, 22, 43 

Optical observatory 

Band 

KVA, Spain 

R 

Xinglong, China 

I 


with a cooled receiver installed at prime focus, with two sym- 
metric off-axis corrugated horn feeds sensitive to left circular 
polarization. The telescope and receiver combination produces 
a pair of approximately Gaussian beams (157 arcsec FWHM), 
separated in azimuth by 12.95 arcmin. The receiver has a cen- 
ter frequency of 15.0 GHz, a 3.0 GHz bandwidth, and a noise- 
equivalent reception bandwidth of 2.5 GHz. Measurements are 
made using a Dicke-switched dual-beam system, with a second 
level of switching in azimuth to alternate source and sky in each 
of the two horns. Calibration is referred to 3C286, for which a 
flux density of 3.44 Jy at 15 GHz is assumed (Baars et al. |l977| l. 
Details of the observations, calibration, and analysis are given 
by ( Richards et al,|20TT| ). 


3.1.7. RATAN 


A six-frequency broadband radio spectrum was obtained with 
the RATAN-600 radio telescope in transit mode by observing si- 


multaneously at 1.1, 2.3, 4.8, 7.7, 11.2, and 21.7 GHz (Parijskij 
1993; [Berlin & Friedmanp996|>. Data were reduced using the 


RATAN standard software FADPS (Flexible Astronomical Data 
Processing System) reduction package ( jVerkhodanov 1997 1 . 
The flux density measurement procedure at RATAN-600 is de- 
scribed by ( jAliakberov et al. 1 1985[ ). 


3.1.9. V LA 

The Very Large Array (VLA) and, since Spring 2010, the 
Expanded VLA (EVLA), observed a subset of the sources as 
simultaneously as possible. Most of the VLA and EVLA runs 
were 1-2 hour chunks of time. In a one-hour chunk of time, in 
addition to flux calibrators and phase calibrators, typically 5-8 
Planck sources were observed. In many cases, VLA flux-density 
and phase calibrators were themselves of interest, since they 
were bright enough to be detected by Planck. For these bright 
sources, the integration times could be extremely short. 

Integration times were about 30 s at 4.86 GHz and 8.46 GHz, 
100 s at 22.46 GHz, and 120 s at 43.34 GHz. All VLA/EVLA 
flux density measurements were calibrated using standard val- 
ues for one or both of the primary calibrator sources used by 
NRAO, 3C48 or 3C286, and the u-v data were flagged, cal- 
ibrated, and imaged using standard NRAO software (AIPS or 
CASA). It is important to bear in mind that the VLA and EVLA 
were in different configurations at different times in the several 
months duration of the observations. As a consequence, the an- 
gular resolution changed. In addition, for a given configuration, 
the resolution was much finer at higher frequencies. For that rea- 
son, sources that showed signs of resolution in any VLA config- 
uration or at any VLA frequency have been carefully flagged. 


3.1.10. KVA 


3.1.8. UMRAO 


Centimetre-band observations were obtained with the University 
of Michigan 26-m prime focus paraboloid equipped with ra- 
diometers operating at central frequencies of 4.8, 8.0, and 
14.5 GHz. Observations at all three frequencies utilized rotating 
polarimeter systems permitting both total flux density and linear 
polarization to be measured. A typical measurement consisted 
of 8 to 16 individual measurements over a 20-40 minute time 
period. Frequent drift scans were made across stronger sources 
to verify the telescope pointing correction curves; and observa- 
tions of program sources were intermixed with observations of 
a grid of calibrator sources to correct for temporal changes in 
the antenna aperture efficiency. The flux scale was based on ob- 
servations of Cassiopeia A (Baars et al. [1977). Details of the 
calibration and analysis techniques are described by |Aller et al. 
( fl985| l. 


Optical observations were made with the 35 cm KVA (Kunliga 
Vetenskapsakademiet) telescope at La Palma, Canary islands. 
All observations were made through the R-band filter (/l e ff = 
640 nm) using a Santa Barbara ST-8 CCD camera with a gain 
factor of 2.3 e ~ / ADU and readout noise of 14 electrons. Pixels 
were binned 2x2 pixels giving a plate scale of 0.98 arcsec/pixel. 
We obtained 3-6 exposures of 180 s per target. The images were 
reduced in the standard way of subtracting the bias and dark 
frames and dividing by twilight flat-fields. The fluxes of the tar- 
get and 3-5 stars in the target field were measured with aper- 
ture photometry and the magnitude difference between the target 
and a primary reference star in the same field was determined. 
Using differential mode makes the observations insensitive to 
variations in atmospheric transparency and accurate measure- 
ments can be obtained even in partially cloudy conditions. The 
R-band magnitude of the primary reference star was determined 
from observations made on photometric nights, using compari- 
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son stars in known blazar fields as calibrators (Fiorucci & Tosti 


1996| Fiorucci et al.|1998[|Raiteri et al.|1998|[Villata et al.|1 998 1 


Nilsson et al.|2007 1 and taking into account the color term of the 
R-band filter employed. After the R-band magnitude of the pri- 
mary reference star was determined, the object magnitudes were 
computed from the magnitude differences. At this phase we as- 
sumed V - R = 0.5 for the targets. Several stars in the field were 
used to check the quality of the photometry and stability of the 
primary reference. The error bars of the magnitudes include the 
contribution from both measurement and calibration errors. 


3.1 .1 1 . Xinglong 

The monitoring at Xinglong Station, National Astronomical 
Observatories of China, was performed with a 60/90 cm f/3 
Schmidt telescope. The telescope is equipped with a 4096x4096 
E2V CCD, which has a pixel size of 1 2 pm and a spatial res- 
olution of l'/3 pixel -1 . The observations were made with an 
intermediate-band filter, the I filter. Its central wavelength and 
passband width are 6685 A and 5 14 A, respectively. The expo- 
sure times were mostly 120 s but ranged from 60 to 180 s, de- 
pending on weather and lunar phase. 


3.2. Planck microwave data 

Planck ( |Tauber et al.|2010| |Planck Collaboration||201 la| is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensi- 
tivity and angular resolution from 31' to 5'. Full sky cover- 
age is attained in ~ 7 months. The Low Frequency Instrument 


(LFI; Mandolesi et al.|2010[ |Bersanelli et al.|2010[ , |Mennella et 
aT7||201 1[ > covers the 30, 44, and 70 GHz bands with amplifiers 
cooled to 20 K. The High Frequency Instrument (HFI; Lamarre 
et al. 2010| |Planck HFI Core Team|201 laj > covers the 100, 143, 
217, 353, 545, and 857 GHz bands with bolometers cooled to 
0.1 K. Polarization is measured in all but the highest two bands 
( Leahy et al.|2010 Rosset et al,|2010[ i. A combination of radia- 
tive cooling and three mechanical coolers produces the temper- 
atures needed for the detectors and optics (Planck Collaboration 


2011b 


brate t 


2011b 


Two Data Processing Centers (DPCs) check and cali 
re data and make maps of the sky (Planck HFI Core Team 
Zacchei e t ak]|20lT| . Planck ' s sensitivity, angular reso- 


lution, and frequency coverage make it a powerful instrument 
for Galactic and extragalactic astrophysics as well as cosmol- 
ogy. Early astrophysics results are given in Planck Collaboration, 
201 lh-z. 

The Early Release Compact Source Catalog (ERCSC, 
|Planck Collaboration|201 lc] > contains all sources, both Galactic 
and extragalactic, detected with high confidence over the full 
sky during the period between August 12, 2009 and June 6, 
2010 (corresponding to Planck operational days 91 to 389). 
The ERCSC only contains average intensity information for the 
sources. However, many of the sources were observed more than 
once during the time period spanned by the ERCSC, and some 
of the Swift observations used for this paper have been carried 
between June and October 2010. Therefore, in order to have 
simultaneous data, we have produced separated maps for the 
first (OD 91-274), the second (OD 275^-56), and the begin- 
ning of the third Planck survey (OD 457-550) through the LFI 
and HFI pipelines described in Zacchei et al. ( |201 1 1 and Planck 
|HFI Core Team| ( |201 lb[ >, and we have extracted the flux den- 
sities from each map using IFACMEX, an implementation of 


the Mexican Hat Wavelet 2 (MHW2) algorithm available at the 
LFI DPC. The MHW2 tool has been extensively used to detect 
point-like objects in astronomical images, both with simulations 
from various experiments and data from the WMAP, Planck, and 
Herschel satellites ( Gonzalez-Nuevo et al.|2006 ; Lopez-Caniego 


|et al.|2006l |2007[ |Massardi et al.|2009| l. This wavelet is defined 
as the fourth derivative of the two-dimensional Gaussian func- 
tion, where the scale of the filter R is optimized to look for the 
maximum in the signal-to-noise ratio of the sources in the fil- 
tered map. In practice, the IFCAMEX code, our implementation 
of the MHW algorithm, deals with FITS maps in Healpix for- 
mat ( jGorski et al.|2005) and can be used to detect sources in the 
whole sky or at the position of known objects. For this analysis 
we have looked for objects above the 4cr level at the positions 
corresponding to the 105 sources of our sample. The results of 
the analysis of Planck simultaneous data are reported in Table [6] 
where columns 1 and 2 give the source name, columns 3 and 4 
give the observation start and end times, and columns 5-13 give 
the flux densities in units of Jy at 30, 44, 70, 100, 143, 217, 353, 
545, and 857 GHz. 

Owing to source variability, we do not expect these simulta- 
neous flux densities to be the same as the time-averaged ERCSC 
measurements, except in the case of the sources that have been 
observed only once during the ERCSC time range and for which 
we have estimated the Planck flux densities, simultaneous with 
the Swift observation, in the same period. We have verified that, 
for the sources fulfilling these requirements, the flux densities 
extracted for this paper are in good agreement with those of the 
ERCSC. 


3.3. Swift optical, UV and X-ray data 


The Swift Gamma-Ray-Burst (GRB) Explorer (Gehrels et al. 
[2004b| > is a multi-frequency space observatory devoted to the 
discovery and rapid follow up of GRBs. There are three instru- 
ments on board the spacecraft: the UV and Optical Telescope 
(UVOT, Roming e t al~| |2005) , the X-Ray Telescope (XRT, 
|Burrows et al.|2005| ) sensitive in the 0.3-10.0 keV band, and the 
Burst Alert Telescope (BAT, Barthelmy et al.|2005 i sensitive in 
the 15-150keV band. Although the primary scientific goal of 
the satellite is the observation of GRBs, the wide frequency cov- 
erage is suitable for blazar studies, because it covers the region 
where the transition between synchrotron and inverse-Compton 
emission usually occurs. 

When not engaged in GRB observations. Swift is available 
for Target of Opportunity (ToO) requests, and the Swift team de- 
cided to devote an average of three ToO observations per week 
to this project for simultaneous observation of blazars. 


3.3.1. UVOT 

The Swift UVOT telescope can produce images in each of its 
six lenticular filters (V, B, U, UVW1, UVM2, and UVW2). 
However, in an effort to reduce the use of mechanical parts after 
several years of orbital operations, observations are carried out 
using only one filter, unless specifically requested by the user. 
Thus images are not always available for all filters. 

The photometry analysis of all our sources was performed 
using the standard UVOT software distributed within the 
HEAsoft 6.8.0 package and the calibration included in the lat- 
est release of the “Calibration Database”. A specific procedure 
has been developed at the ASI Science Data Center (ASDC) 
to process all ToO observations requested for the blazar sam- 
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pie. Counts were extracted from aperture of 5" radius for all 
filters and converted to fluxes using the standard zero points 
( |Poole et al,||2008) . The fluxes were then de-reddened us- 
ing the appropriate values of E(B - V ) for each source taken 
from ( Schlegel et al.|1998] > with AxlE{B - V) ratios calculated 
for UVOT filters using the mean Galactic interstellar extinc- 
tion curve from ( |Fitzpatrick| 1999) ). No variability was detected 
within single exposures in any filter. The processing results were 
carefully validated including checks for possible contamination 
from nearby objects within the source and background apertures. 
Some sources, such as 3C 273 and NGC 1275, needed special 
analysis, and results for some other sources had to be discarded. 

The results of UVOT data analysis are summarized in Table[7] 
where columns 1 and 2 give the source name, columns 3 and 4 
give the observation date and the Swift observation ID, and the 
remaining columns give the magnitudes in the six UVOT filters 
with errors. 


3.3.2. XRT 


The Swift XRT is usually operated in “auto-state” mode which 
automatically adjusts the CCD read-out mode to the source 
brightness, in an attempt to avoid pile-up ( jBurrows et al. 
2005 1 Hill et ah) 2004 1. As a consequence, some of the data 
were collected using the most sensitive Photon Counting (PC) 
mode while Windowed Timing (WT) mode was used for bright 


The XRT data were processed with the XRTDAS software 
package (v. 2.5.1, Capalbi et al. (j2005j)) developed at the ASI 


Science Data Center (ASDC) and distributed by the NASA High 
Energy Astrophysics Archive Research Center (HEASARC) 
within the HEASoft package (v. 6.9). Event files were calibrated 
and cleaned with standard filtering criteria using the xrtpipeline 
task and the latest calibration files available in the Swift CALDB. 
Events in the energy range 0.3-10 keV with grades 0-12 (PC 
mode) and 0-2 (WT mode) were used for the analysis. 

Events for the spectral analysis were selected within a cir- 
cle of 20 pixel (~ 47") radius, which encloses about 90% of the 
PSF at 1 .5 keV ( jMoretti et al,|2005) , centered on the source po- 
sition. For PC mode data, when the source count rate is above 
~ 0.5 counts/s data are significantly affected by pile-up in the in- 
ner part of the Point Spread Function (PSF). For such cases, after 
comparing the observed PSF profile with the analytical model 
derived by (Moretti et al. 2005) 1, we removed pile-up effects by 
excluding events detected within up to 6 pixels from the source 
position, and used an outer radius of 30 pixels. The value of the 
inner radius was evaluated individually for each observation af- 
fected by pile-up, depending on the observed source count rate. 

Ancillary response files were generated with the xrtmkarf 
task applying corrections for the PSF losses and CCD defects. 
Source spectra were binned to ensure a minimum of 20 counts 
per bin to utilize they 2 minimization fitting technique. We fitted 
the spectra adopting an absorbed power-law model with pho- 
ton index When deviations from a single power law model 
were found, we adopted_a log-parabolic law of the form F(E) = 
KE (a+bu,gE) (jMassaro et alj|2004j» which has been shown to fit 


the X-ray spectrum of blazars of the HSP type well (e.g. Giommi 
et al.||2005) Tramacere et al.|2009) . This spectral model is de- 
scribed by only two parameters: a, the photon index at 1 keV, 
and b, the curvature of the parabola. For both models the amount 
of hydrogen-equivalent column density (Ah) was fixed to the 
Galactic value along the line of sight ( Kalberla et al.|2005 1 . For 
a fraction of the sources fitted with a power-law model (~15%) 


we found evidence of an absorption excess at low energies and 
the hydrogen column density Ah parameter was left free. 

The results of the spectral fits with a power-law model and 
Galactic Ah are shown in Table[8]where columns 1 and 2 give the 
source name, column 3 gives the Swift observation date, column 

4 gives the best fit photon index T x , column 5 gives the Galactic 
Ah, columns 6 and 7 give the 0. 1-2.4 and 2-10 keV X-ray fluxes, 
column 8 gives the value of the reduced y 2 , and column 9 gives 
the number of degrees of freedom. 

In Table [9] we report data obtained using a log-parabola to 
describe the spectrum model. Column 1 and 2 give the source 
name, column 3 gives the Swift observation date, column 4 and 

5 give the log parabola parameters a and b, column 6 gives the 
Galactic Ah, columns 7 and 8 give the 0. 1-2.4 and 2-10 keV 
X-ray fluxes, column 9 gives the value of the reduced y 2 - and 
column 10 gives the number of degrees of freedom. 


Finally, in Table 10 we report the results of the spectral fits 
that have been performed leaving the hydrogen column den- 
sity Ah as a free parameter. Columns 1 and 2 give the source 
name, column 3 gives the Swift observation date, column 4 gives 
the best fit photon index T x , column 5 gives the estimated Ah, 
columns 6 and 7 give the 0. 1-2.4 and 2-10 keV X-ray fluxes, 
column 8 gives the value of the reduced y 2 , and column 9 gives 
the number of degrees of freedom. 


3.4. Fermi-MT y-ray data 

The Large Area Telescope (LAT) on-board Fermi is an electron- 
positron pair conversion telescope sensitive to y-rays of ener- 
gies from 20MeV to > 300 GeV. The Fermi-LAT consists of 
a high-resolution silicon microstrip tracker, a Csl hodoscopic 
electromagnetic calorimeter and an anticoincidence detector for 
charged particle background identification. A full description of 
the instrument and its performance can be found in |Atwood et al.| 
(2009 1 . The large field of view (~2.4sr) allows the LAT to ob- 
serve the full sky in survey mode every 3 hr. The LAT point 
spread function (PSF) depends strongly on both the energy and 
the conversion point in the tracker, but less on the incidence an- 
gle. 

The LAT y-ray spectra of all AGN sources are studied in 
|Abdo et al.[ ( |2010b| l based on 1 1 months of Fermi-LAT data. 
Here we derive the y-ray spectra of the blazars for which we 
build the simultaneous SEDs, integrating for two weeks encom- 
passing the whole duration of the Planck observations. 

The Fermi-LAT data considered for this analysis cover the 
period from August 4, 2008 to November 4, 2010 and have been 
analyzed using the standard Fermi - LAT ScienceTools software 
packag^] (version v9rl6 and selecting for each source only pho- 
tons above 100 MeV, belonging to the diffuse class (Pass6 V3 
IRF; |Atwood et al. 2009) which have the lowest background con- 
tamination. For each source, we selected only photons within a 
15° region of interest (Rol) centered around the source itself. In 
order to avoid background contamination from the bright Earth 
limb, time intervals where the Earth entered the LAT Field of 
View (FoV) were excluded from the data sample. In addition, 
we have excluded observations in which the source under study 
was viewed at zenith angles larger than 105° where Earth’s at- 
mospheric y-rays increase the background contamination. The 
data were analyzed with a binned maximum likelihood tech- 
nique (Mattox et al. 1996 1 using the analysis software gtlike de- 


7 http : // fermi . gsf c . nasa . gov/ssc/data/ analysis/ 
documentation/Cicerone 
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veloped by the LAT tean^j A model accounting for the diffuse 
emission and nearby y-ray sources is included in the fit. 

The diffuse foreground, including Galactic interstellar emis- 
sion, extragalactic y-ray emission, and residual CR background, 
has been modelled with gll_iem_v 02^] for the Galactic diffuse 
emission and isotropic iern v()2 for the extragalactic isotropic 
emission. Each source under study has been modelled with a 
power law function 


dN 

dE 


N(T+l)E r 
F r+i _ F r+i 

^max ^ m i n 


( 1 ) 


where both the normalization factor N and the spectral index T 
are left free in the model fit. The model also includes all the 
sources within a 20° Rol included in Fermi-LAT one year cat- 
alog Abdo et al. ( 20 1 0c) > and modelled using power law func- 
tions. If a source included in the model is a pulsar belonging 
to the Fermi-LAI pulsar catalog Abdo et al. ( 2010dj >, we have 
modelled the spectrum with a power law with exponential cut- 
off using the spectral parameters in the pulsar catalog. 

For the evaluation of the y-ray SEDs, the whole energy range 
from 100 Me V to 300 GeV is divided into two equal logarithmi- 
cally spaced bins per decade. In each energy bin the standard gt- 
like binned analysis has been applied assuming power-law spec- 
tra with photon index = -2.0 for all the point sources in the 
model. Assuming that in each energy bin the spectral shape can 
be approximated by a power law, the flux of the source in all se- 
lected energy bins is evaluated, requiring in each energy bin a TS 
greater than 10. If the TS is lower than 10, an upper limit (UL) 
is evaluated in that energy bin. Only statistical errors for the fit 
parameters are shown in the plots. Systematic errors, due mainly 
to uncertainties in the LAT effective area derived from the on- 
orbit estimations, are < 5% near 1 GeV, 10% below 0.1 GeV, 
and 20% above 10 GeV. 

For each source, we have considered the three different inte- 
gration periods for the y-ray data: 


- Simultaneous observations : data accumulated during the pe- 
riod of Planck observation of the source. As the Planck in- 
struments point in slightly different directions and the field of 
view depends on the frequency of observation, a typical ob- 
servation covering all Planck channels takes about one week, 
the exact integration time depending on the position of the 
source. 

- Quasi-simultaneous observations: data integrated over a pe- 
riod of two months centered on the Planck observing period 
of the source. 

- 27 month Fermi-LAT integration: data integrated over a pe- 
riod of 27 months from August 4, 2008 to November 4, 2010, 
i.e., the entire Fermi-LAT data set available for this paper. 

Tables m and [12] give a summary of the y-ray detections 
(TS > 25) in all our samples. The fraction of sources detected by 
Fermi- LAT during the simultaneous integrations is not very large 
and varies from ~ 40% in the Fermi- LAT sample to just ~ 20% 
in the soft X-ray selected sample. We note that even considering 
all the Fermi-LAI data available at the time of writing (27 month 
integration), a sizable fraction of the blazars in the radio, soft and 
hard X-ray selected samples have not been detected. 


8 http : //f ermi . gsf c . nasa . gov/ ssc/data/analysis/ 
documentation/Cicerone/Cicerone_Likelihood 
v http : //f ermi . gsf c . nasa . gov/ ssc/data/access/lat/ 
BackgroundModels . html 


Table 11. Summary of y-ray detections with significance TS > 
25. 


Sample 

No. of detected sources 
simult. 2 months 27 months 

Sources 
in sample 

Fermi - LAT 

18(36%) 

40 (80%) 

50 (100%) 

50 

Swift - BAT 

9 (26%) 

12 (35%) 

27 (79%) 

34 

ROSAT/R ASS 

10 (23%) 

15 (35%) 

24 (56%) 

43 

Radio 

22 (21%) 

38 (37%) 

78 (75%) 

104 


Table 12. Statistics of y-ray detections (TS > 25) in the 27 
month Fermi-LAT data set. 


Sample 

No. of detected sources 
FSRQs BL Lacs Uncertain 

Fermi-LAT 

28 (100%) 

14(100%) 

8 (100%) 

Swift-BAI 

17 (63%) 

7(100%) 

3 (50%) 

ROSAT/RASS 

8 (53%) 

14 (88%) 

2 (17%) 

Radio 

48 (72%) 

16(100%) 

9 (64%) 


Detailed results of the Fermi-LAT analysis are given in 
Tables 13]- 18 where the observed fluxes or upper limits are given 
in six or three energy bands depending on the source brightness. 

Two sources (PKS 0548-322 and NGC 7213) appear as sig- 
nificant y-ray detections in our 27 month data set although they 
have not been included in any of the Ferm i-LAT catalogues pub- 


lished so far ( Abdo et al.||2009b 2010c andp}. These should 


therefore be considered as new y-ray detections. 


4. The importance of simultaneity 


Blazars are, by definition, highly variable sources. It is there- 
fore important to use simultaneous multi-frequency data to build 
SEDs for comparison with theoretical models. In this section we 
compare our measurements with data taken from the literature in 
order to derive an estimate of the uncertainties introduced by the 
use of non-simultaneous data in different parts of the spectrum. 

Figure [T| plots the Planck flux density at 44 GHz presented 
in this paper versus the WMAP flux density at 41 GHz from the 
WMAP point source catalogs ( Bennett et al.|2003] Wright et al. 
2009a). Some scatter is present, but most of the points lie be- 
tween the two solid lines indicating a factor of two variability. 

Figure [2] plots the X-ray fluxes of the sources observed by 
Swift simultaneously with Planck (see Table[8]l against the X-ray 
fluxes of the same sources from the BZCAT catalog (Massaro 
et al.1 |2009[ |20 10 1 . In this case a large scatter is present, with 
variations of over a factor of ten. 

Figure [3] shows the Fermi-LAT y-ray fluxes of our sources 
measured simultaneously with Planck plotted against their y-ray 
fluxes in the Fermi-LAI 1FGL catalog ( |Abdo et al.|2010c i. As 
in the X-ray sample, a scatter with variations larger than a factor 
of ten is observed. 

We conclude that SEDs built with non-simultaneous data 
suffer from uncertainties in the microwave region that are rel- 
atively modest and generally limited to about a factor of two, 
while the high energy part of the spectrum (X-ray and y-ray) 
is much more affected, with uncertainties due to flux variations 


10 http : //heasarc . gsfc . nasa . gov/W3Browse/ all/ 
fermilpsc.html 
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WMAP flux density (41 GHz, mjy) 

Fig. 1. The Planck 44 GHz flux density of the sources in our sam- 
ple is plotted against the 41 GHz flux density from the WMAP 
5-year catalog (81 sources). The three solid lines represent equal 
flux densities (i.e., no variation) and a factor of two variability 
above or below the equal flux level. Almost all the points lie 
between the factor of two variability lines. 



Fig. 3. The Fermi- LAT y- ray flux of the sources in our samples 
detected during the simultaneous integration with Planck is plot- 
ted against the flux reported in the Fermi -LAT 1-year catalog. 
The three solid lines represent equal fluxes (i.e., no variation) 
and a factor of two variability above or below the equal flux 
level. Note that several points are outside the factor of two vari- 
ability lines, revealing variability of up to about a factor ten. 



Simultaneous X-ray flux (0.1-2.4KeV, erg cnv 2 s~') 

Fig. 2. The Swift X-ray (0. 1-2.4 keV) flux of the sources in 
our sample measured simultaneously with Planck is plotted 
against the 0. 1-2.4 keV flux reported in the BZCAT catalog (83 
sources). The three solid lines represent equal fluxes (i.e., no vari- 
ation) and a factor of two variability above or below the equal 
flux level. Note that several points are outside the factor of two 
variability lines, revealing variability of up to about a factor ten. 


of up to a factor of ten or more. Of course the same uncertain- 
ties apply when searching for correlations in non-simultaneous 
multi-frequency data. 


5. Spectral Energy Distributions 

We have constructed the SEDs of all the blazars in our samples 
from the simultaneous multi-frequency data described above us- 


ing the ASDC SED Builder, an on-line service developed at the 
ASI Science Data Center (ASDC^ 11 (Strattaetal. 2011 1 . This is a 


WEB-based tool that allows users to build multi-frequency SEDs 
combining data from local catalogs and external services (e.g., 
NED, SDSS, USNO) with the user’s own data. The tool con- 
verts observed fluxes or magnitudes into de-reddened fluxes at a 
given frequency using standard recipes that take into account the 
instrument response and assumed average spectral slopes. The 
SED builder can display SEDs both in flux and in luminosity (if 
redshift information is given); it also provides useful tasks such 
as the overlay of templates for blazar host galaxies and nuclear 
optical emission (blue-bump), and allows users to compare the 
SED with models including one or more SSC components. 

The SEDs of all the sources in our samples are shown in 
Figs. 25 -42] Red points represent strictly simultaneous multi- 
frequency data; green points represent y-ray data integrated over 
a period of two months centered on the Swift/Planck observa- 
tions, or ground-based data taken quasi-simultaneously; blue 
points represent y-ray data integrated over the full period of 27 
months. In the few cases where no Swift simultaneous observa- 
tions could be obtained we plot only Planck, Fermi-LAT, and 
ground-based data. Two-sigma upper limits are indicated by ar- 
rows. 


5.1. Separating the non-thermal/jet-related radiation from 
QSO accretion and host galaxy emission 

In the following we use the simultaneous SEDs of Figs |25]j42l 
to determine some parameters that can constrain the physical 
mechanisms powering blazars. However, before doing so we 
must identify and separate radiation that is unrelated to the non- 
thermal, relativistically amplified emission from the jet of the 
blazars; that is, radiation from accretion onto the central black 


11 http://tools.asdc.asi.it/SED/ 
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Log (frequency), Hz 


Fig. 4. Top panel: The SDSS template of Vanden Berk et al. 


2001 


for the broad-line and thermal emission from a QSO. Bottom 
panel : The giant elliptical galaxy template of Mannucci et al. 
2001 See text for details. 


hole and from the host galaxy (see, e.g., Perlman et al.||2008) . 
To do so we estimated the contamination from the host galaxy 
assuming that all blazars are hosted by giant elliptical galaxies 


(e.g., Kotilainen et al.|[l998| Nilsson et al.||2003j Leon-Tavares 


201 1[) with absolute magnitude of Mr = -23.7. As Scarpa 


et al. 


et al. 


( 2000 1 and Urry et al. (|2000) have demonstrated, this value 


is within 1 magnitude of the observed values in a sample of 
over 100 BL Lacs observed with the Hubble Space Telescope. 
For the spectral shape we used the elliptical galaxy template of 
Mannucci et ak][2001| (bottom panel of Fig. [4), which is based 


on low resolution observations of a number of nearby galaxies 
in the wavelength range 0.12 -2.4 n m and is a good match to 
the predictions of spectrophotometric models for giant ellipti- 
cals (Mannucci et ai .[200T l 

The radiation due to accretion was estimated from the com- 
posite optical spectrum built by Vanden Berk et al. ( |2001| > using 


over 2200 optical spectra of radio-quiet QSOs taken from the 
SDSS database ( |York et aTp OOO ) (top panel of Figj4|), and the 
expected soft X-ray emission of radio quiet AGN from |Grupe 
let al.l2010t 

The ratio of optical to soft X-ray light has been known to be 
a function of optical luminosity since the time of the Einstein ob- 
servatory (Avni & Tananbau m|1986j l. Recently this dependence 
has been confirmed using simultaneous optical and soft X-ray 
data from Swift (|Grupe et al.|2010|l and XMM-Newton ( Vagnetti 


et al. 


et al. 


2010[>. In this paper we use the relationship given by Grupe 

( 2010 ): 


«UV-X(Radio-quietQSO) = 0. 1 14 log(L 2500 - 1.177 


( 2 ) 


where L 05(|(| A is the rest-frame luminosity of the thermal emis- 
sion at 2500 A in units of W Hz -1 and auv-X(radio-quietQSO) is the 
usual slope between the UV (2500 A) and the soft X-ray (2 keV) 
flux in radio quiet QSOs (e.g., [Vagnetti et a l.|2010 >. 

Examples of the emission from these components unrelated 
to the jet are shown in Fig. [5] which shows the SEDs of Mkn501 
and of Mkn 421, where the optical light is dominated by the host 
galaxy. 


Figs. [ 6 ] and [ 7 ] show the SEDs of BL Lacertae, where part 
of the UV light is thought to be due to disk emission (Raiteri 


et~~aL||2009) , the nearby FSRQ 3C273 (z = 0.158), and the 


high-redshift FSRQs 4C 38.41 (z = 1.814) and lJy 0537-286 
( z = 3.104) where the optical/UV light is heavily, or completely, 
contaminated by radiation coming from the accretion onto the 
central black hole. 

We assess the amplitude of the optical thermal emis- 
sion compared to non-thermal radiation using the parameter 
Q'R-O(Thermai), defined as the spectral slope between the 5 GHz 
radio flux density and the 5000 A optical flux density that can be 
attributed to the blue-bump/disk/thermal emission. This quan- 
tity depends both on the relativistic amplification factor and the 
intrinsic ratio between non- thermal/jet radiation and disk emis- 
sion. In the case of FSRQs the optical spectrum shows emission 
lines by definition, therefore we can constrain «'r -o(Thermai) by 
adjusting the optical thermal emission at the same level as the 
data. When the thermal blue bump is seen directly in the op- 
tical/UV part of the SED (e.g.. Figs. [ 6 ] and [ 7 ) we fit the tem- 
plate of Vanden Berk et al. ( |2001) to the observed data; when 
the disk emission is not obviously visible we set the intensity 
of the blue-bump template just below the observed optical/UV 
emission where the broad emission lines have been detected. 

BL Lacs do not show emission lines in their optical spec- 
trum, so for this class of objects we can only set upper limits on 
(ZR_ 0 (Thermai) ■ We do that by assuming that the template for the 
optical thermal emission is at least one order of magnitude below 
the observed data, that is sufficiently low to hamper any broad 
line detection. The estimation of u-r o(Thermai) relies strongly on 
the quality of the optical data available, in particular on the si- 
multaneous UVOT data. Therefore we define an “optical data 
quality” flag as follows: 


- 0: no simultaneous data available; 

- 1: poor quality (e.g., only one UVOT filter available); 

- 2: good quality (e.g., two or three UVOT filters available); 

- 3: excellent quality (all UVOT filters available). 

In Figure [8] (upper panel) we show the distribution of 
ffR-O(Thermai) for all the FSRQs with excellent optical data for the 
whole sample and for each of our samples independently. The 
distribution of the otr .o(Thermai) upper limits for BL Lac objects is 
shown in the bottom panel. The results suggest a possible differ- 
ence in the o-r o(Therman distribution in different samples. In fact, 
as reported in Table |T9) when only sources with excellent optical 
data are considered we obtain (aR-ocrhermai)) = -0.64 + 0.05 for 
the ROSAT/ RASS sample, to be compared with (o-r oiTherman) = 
-0.73 ± 0.02 for the Fermi- LAT sample. 

In order to assess the possible presence of a thermal com- 
ponent in the X-ray emission, for each source we compare the 
predicted thermal emission from accretion with the actual X-ray 
spectrum, and we also include the uncertainties in the parame- 
ters of Eq.[2 in order to build 1, 2, and 3<x bands. We therefore 
define the fo lowing X-ray thermal contamination flag: 

- 0: no X-ray data available; 

- 1: X-ray emission mostly or entirely due to accre- 
tion/reflection (data agree with the expectations for accretion 
emission within lcr); 

- 2: X-ray emission probably contaminated by the accretion 
component (data agree with with the expectations for accre- 
tion emission within 2cr); 

- 3: X-ray emission mostly of non-thermal origin (data are at 
2-3cr from the expectations for radio quiet QSOs); 
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Fig.5. The SEDs of Mkn501 (left) and Mkn421 (right) showing the expected emission from the host galaxy (giant elliptical) as an 
orange line. The green lines are the best fit to the simultaneous non-thermal data using a third degree polynomial function. See text 
for details. 




Fig. 6. Left : the SED of BL Lacertae showing the expected emission from the host galaxy (just below the observed non-thermal 
radiation, giant elliptical, orange line) and of the blue bump emission (blue line, see also Raiteri et al. 2009 j >. Right : the SED of 


3C273 showing the thermal emission from the blue bump and the expected X-ray emission from accretion (purple and blue lines, 
see text for details). The vertical parallel lines represent the optical window (4000-10000 A). The green lines are the best fit to the 
simultaneous non-thermal data using a third degree polynomial function. 


Table 19. Contamination of the X-ray emission from a thermal component and (ctr. o(Thermai)) estimated for FSRQs in the whole 
sample and for each sample independently, for various thresholds on the quality of the optical data (see the text for a description of 
the flag). 


ALL Fermi - LAT Swift- BAT ROSAT/RASS 


Opt Data 
quality 

cont/ 

total 

% 

(UR-O(Thermal) ) 

cont/ 

total 

% 

(ttR-O(Thermal)} 

cont/ 

total 

% 

(tfR-OlThermal)} 

cont/ 

total 

% 

(«R-0(Thermal)) 

ALL 

10/47 

21 

-0.69±0.02 

3/25 

12 

-0.72±0.02 

2/21 

10 

-0.70±0.03 

6/13 

46 

-0.65±0.04 

>2 

10/44 

23 

-0.69±0.02 

3/22 

14 

-0.73±0.02 

2/20 

10 

-0.70±0.03 

6/13 

46 

-0.65±0.04 

3 

9/40 

23 

-0.69±0.02 

3/22 

14 

-0.73±0.02 

1/16 

6 

-0.69±0.04 

6/12 

50 

-0.64±0.05 


- 4: X-ray emission certainly of non-thermal origin (data are 
more than 3cr away from the expectations for radio quiet 
QSOs). 

Results for all the FSRQs with good or excellent optical data 
are shown in Figure[9] and summarized in Table[l9]for the whole 


sample and for each sample independently. There is a large dif- 
ference between the Fermi - LAT and Swift -BAT samples, where 
< 15% of the sources show a thermal component in the X-ray 
emission, and the ROSAT/RASS sample, where ~ 50% of the 
sources are contaminated. This demonstrates that the thermal 
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Fig. 7. The SEDs of 4C 38.41 (z. = 1.814; left ) and of lJy0537-286 (z = 3.104; right). Simultaneous data are shown in red; 
non simultaneous literature or archival data are shown in light gray. Note that the UVOT data of these medium and high redshift 
objects matches quite well the QSO template of Vanden Berk et al. ( |2001) > that we use to estimate the thermal emission from the 
blue bump ( blue line). The emission from the host galaxy ( orange line ) is very low compared to other components. The observed 
X-ray emission from these sources is more than 3<x above the expected emission from accretion. In both cases the optical light is 
dominated by accretion while the X-rays are from the non-thermal component. 


Radio sample 
| Fermi sample 

BAT sample 
RASS sample 



RO(Thermal) 


Fig. 8. Distribution of aR-o(Thermai) for all the sources and for 
each blazar sample. Only SEDs with good or excellent optical 
data (flag = 2 or 3 , see text) were used. Black solid histograms: 
radio sample; red dashed histograms: Fermi sample; green dot- 
dashed histograms: Swift BAT sample; blue dotted histograms: 
RASS sample. Distributions for BL Lac objects refer to upper 
limits only. 


component in the X-ray emission of blazars cannot be neglected, 
even in bright sources. 


5.2. SED parameter estimation 


We have used the SEDs of all the objects in our samples to esti- 
mate the values of important physical parameters such as v s 


peak’ 


v s ,F(v s ,), v IC , and y ic ,F(v ]c ,,) (see Table 20 1 taking into 

peak v peak 7 peak neak v neak 7 x ° 


peak v peak 


account only the non-thermal radiation and fitting third-degree 
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Fig. 9. Distribution of the thermal emission flag for all the 
FSRQs with good or excellent optical data flag (see text for 
details). Black solid histograms: radio sample; red dashed 
histograms: Fermi- LAT sample; green dot-dashed histograms: 
Swift BAT sample, blue dotted histograms: RASS sample. As all 
BL Lacs showed no contamination of their X-ray emission, they 
are not plotted here. 


polynomials as described in Abdo et al. ( 2010a| > (see Figs. [5] [6] 
and[7]for examples). 

For the sources that were not detected by Fermi- LAT even 
in the 27 month integration we estimated limits to vj^ ak and 

v pe-ik^ y pe-ik) by constraining the polynomial in the high-energy 
part of the SED with the 27 month Fermi- LAT upper limits as 


shown in Fig. 10 


In some HSP BL Lacs with particularly high y s _,, values 


(see, e.g.. Figs. mmm and 


peak 

the Fermi- LAT data alone 
are not sufficient to ensure a good measure of y r ^ ak , as the spectra 
are still rising at the highest Fermi-LAT energies and no simulta- 
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Table 20. Rest-frame (v* eak ) and (v' aak ) values for different samples and classes. 


Sample 

Class 

Synchrotron peak 
No. of sources (log(v/. ent )) 

Inverse-Compton peak 
No. of sources (log(vp ra *)) 


ALL 

45 

13.34±0.12 

42 

22.19+0.10 

Fermi-L AT 

FSRQs 

24 

13.07±0.07 

23 

22.21+0.14 


BL Lacs 

13 

14.11+0.32 

11 

22.35+0.24 


ALL 

26 

13.65+0.26 

22 

21.66+0.28 

Swift - BAT 

FSRQs 

18 

13.16+0.09 

16 

21.45+0.17 


BL Lacs 

5 

15.75+0.81 

3 

23.41 + 1.52 


ALL 

25 

14.27+0.31 

20 

22.15+0.38 

ROSAT/RASS 

FSRQs 

11 

13.02+0.17 

11 

21.43+0.30 


BL Lacs 

11 

15.79+0.28 

6 

24.24+0.37 


ALL 

94 

13.20+0.06 

77 

21.99+0.10 

Radio 

FSRQs 

64 

13.08+0.05 

49 

21.99+0.12 


BL Lacs 

14 

13.89+0.27 

13 

22.26+0.20 


neous TeV data are available. The v IC , values for these sources 

peak 

should therefore be considered as lower limits. 



Log Frequency v (Hz) 


Fig. 10. The SED of the blazar PKS 0003-066, illustrating the 
estimation of the upper limits on vj^ ak and on vJ^ ak F(v l ^a k ) by 
combining the X-ray data with the 27 month Fermi - LAT upper 
limits. 


6. The spectral slope of blazars in the 
radio-microwave region 

While WMAP results are consistent with a single flat spectral in- 
dex in the relatively narrow frequency range 23-94 GHz (Wright 
|et al.|2009b Gold et al.|201 1| ), the blazar spectrum must steepen 
at frequencies closer to the synchrotron peak. Adding Planck 
data and simultaneous ground-based observations at centimetre 
wavelengths to the WMAP data improves the spectral coverage, 
and allows us to probe the spectral shape of blazars over the 
much wider frequency range ~ 1 GHz to ~ 1 THz. 

We have studied the low frequency (LF) and the high fre- 
quency (HF) regions of the centimetre to sub-millimetre blazar 
spectra separately in order to search for differences in the spec- 
tral index a and to determine the frequency at which the spectral 
index changes. We have fitted the two frequency regions inde- 
pendently with power laws to estimate the spectral indices at low 
frequencies (cflf, for v < VBreak) and at high frequencies (othf. 


for v > v Brea k), assuming a range values of the break frequency, 
VBreak. from 30 to 100 GHz. In Fig.[TT]we show the distributions 
of QfLF and ant for VB rea k = 70 GHz and 100 GHz; the blazar 
spectra steepen from rr LF - 0 to cthf 0.65. 

In order to verify the robustness of the results we have re- 
peated the analysis imposing two minimum numbers of inde- 
pendent frequencies needed to perform the fit, namely 3 and 5, 
and the results are consistent (Table [2T|. 

We have also analysed the different classes of blazars - Flat 
Spectrum Radio Quasars (FSRQ), BF Fac objects (BF Fac), 
and AGN of unknown classification (AGN) - separately, and we 
have found a difference at a level of about 3cr in ant, which is 
~ -0.7 for FSRQ but ~ -0.5 for BF Facs. 

The spectral steepening and the difference in cthf are visible 
in Figure [12] where we show «-|> vs. ant for the sources that 
meet the requirements on the minimum number of independent 
data at both low and high frequency. 

Our results are in complete agreement with the findings of 
Planck Collaboration] ^!) 1 Te for the radio selected sample, i.e., 
flat spectral index at low frequency and a H f 0.6 at high fre- 

quency with a break frequency ~ 70 GHz, suggesting that this 
is a general feature of all blazars regardless of the selection cri- 
teria. There is also general agreement with the results presented 
in lPlanck Collahorationl201 ldi for all the sources in the Planck 
ERCSC catalogs, and the slight discrepancy in the spectral index 
estimated at low frequency can be explained by the fact that, un- 
like Planck Collaboration 201 ld[ we have included the ground- 
based observations at 5-30 GHz. 


7. Searching for correlations between fluxes in 
different energy bands 


In this section we use our simultaneous and average multi- 
frequency data to test for the presence of correlations between 
fluxes measured in different energy bands. As some of the 
blazars in our samples have not been detected by Fermi - FAT 
or by Planck, we estimate the significance of the correlations us- 
ing the ASURV code Rev 1.2 ( Lavalley et al.|1992]>, wh ich takes 
account of flux limits as described in|Isobe et al.|(|1986]. 


7. 1 . Microwave vs. X-ray 

Figure [13] shows the Planck 143 GHz flux density versus the si- 
multaneous Swift-X RT X-ray flux for all the sources where the 
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Fig. 11. Radio/microwave low frequency (v < 70 GHz, LF; left side) and high frequency (v > 70 GHz, HF; right side) spectral 
index distributions of FSRQs and BL Lacs in our samples for the case v Break = 70 GHz. While the distributions of low frequency 

spectral indices are very similar in both types of blazars with (au) 0.1, the distributions of the high frequency slopes indicate 

that FSRQs ((qihf) = -0.73 ± 0.04) might be somewhat steeper than BL Lacs ((cthf) = -0.51 ± 0.07). A KS test gives a probability 
of less than 3% that the HF distributions for FSRQs and BL Lacs come from the same parent distribution. 

Table 21. Radio LF and HF spectral index distributions for different values of the break frequency v Break and of the minimum 
number of frequency bands considered for the fit. 


V Break 

n. of frequencies 

Class 

Low Frequency 
n. of sources (au) 

o"lf 

High Frequency 
n. of sources (cthf) 

a- hf 

30 

3 

ALL 

47 

0.00± 0.04 

0.27 

69 

-0.5 1± 0.04 

0.30 

30 

5 

ALL 

46 

0.00± 0.04 

0.27 

57 

-0.49± 0.03 

0.26 

44 

3 

ALL 

47 

-0.02± 0.04 

0.25 

66 

-0.64± 0.03 

0.26 

44 

5 

ALL 

47 

-0.02+ 0.04 

0.25 

47 

-0.62± 0.03 

0.22 

70 

3 

ALL 

63 

-0.08± 0.03 

0.28 

66 

-0.67± 0.03 

0.28 

70 

5 

ALL 

47 

-0.02± 0.04 

0.25 

35 

-0.65± 0.04 

0.22 

100 

3 

ALL 

65 

-0.1 1± 0.03 

0.27 

48 

-0.56± 0.04 

0.26 

100 

5 

ALL 

47 

-0.02± 0.04 

0.24 

15 

-0.52± 0.05 

0.19 

70 

3 

FSRQ 

39 

-0.1 1± 0.04 

0.28 

40 

-0.73± 0.04 

0.27 

70 

3 

BL Lac 

12 

-0.08± 0.08 

0.29 

14 

-0.51± 0.07 

0.27 

70 

3 

Unc. Type 

12 

-0.02± 0.08 

0.29 

12 

-0.64± 0.07 

0.25 


X-rays are expected to be due to the inverse-Compton compo- 
nent (that is all LSP blazars) and not significantly contaminated 
by X-ray emission that is unrelated to the jet such as that due 
to the accretion process (see Sect. 5.1 for details). Sources that 
were not detected by Planck are p otted as upper limits; all the 
sources were detected in the X-ray band. 

A correlation, although with some scatter, is clearly present. 
The Spearman rank coefficients (p) and the corresponding prob- 
abilities that the observed correlation is the result of chance are 


given in Table 22 Results obtained using Planck flux densities 
at other frequencies from 30 to 217 GHz are similar and are not 
shown here. 


7.2. Microwave vs. y-ray 

The relationship between radio/microwave and y- ray fluxes is a 
topic that has been addressed several times in the literature; gen- 
erally a positive correlation between the radio and y-ray fluxes 
is found, though with a large scatter (e.g., Kovalev et al.|2009~ 


2010; Mahony et al. 2010 Peel et al. 2010 

Linford et al. 2011; 

Ackermann et al. 201 1 ; Leon-Tavares et al. 

2011). However, in 


most cases radio and non-simultaneous y- ray data for sources 
detected in both energy bands are compared; here for the first 
time we present simultaneous microwave and y-ray data and we 
take into account upper limits. 


|Leon-Tavares et al~~||20 1 Oj [Giroletti et al.|2010| IGhirlanda et al.| 


The top panel of Figure 14 shows the simultaneous Ferini- 
LAT y-ray flux plotted versus the Planck flux density at 
143 GHz. Sources with TS < 25 in the Fermi-L AT data and 
below 4cr in Planck maps are shown as upper limits. The 
plot shows a clear trend; but from Table 23 which gives the 
Spearman’s p correlation parameter and the probability P that 
the observed level of correlation is due to chance, we see that 
the level of significance of the correlation is never very high (P 
is of the order of a few percent and never lower than 0.05%), 
especially in the hard X-ray and y-ray selected samples and for 
BL Lacs. This result is partly due to the large number of upper 
limits in the simultaneous y-ray flux, which was estimated typi- 
cally over a period of one week. To improve the statistics, in the 
middle panel of Fig. [14] we show the same plot using the Fermi- 
LAT y-ray flux integrated over a two month period centered on 
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Table 22. Spearman correlation parameters for 143 GHz flux density vs. X-ray flux. 


Class Fenni-LRT sample Swift-B AT sample ROSAT/RASS sample Radio sample 



No. of 
sources 
(a)/(b) 

P 

P{%) 

No. of 
sources 
(a)/(b) 

P 

P(%) 

No. of 
sources 
(a)/(b) 

P 

P(%) 

No. of 
sources 
(a)/(b) 

P 

P{%) 

All 

35/3 

0.72 

<0.01 

18/4 

0.50 

2.34 

9/10 

0.70 

0.30 

49/3 

0.51 

0.03 

FSRQs 

20/2 

0.75 

0.06 

16/3 

0.58 

1.40 

5/3 

0.89 

1.85 

33/3 

0.57 

0.08 

BL Lacs 

8/1 

0.95 

0.72 

2/0 

...* 


1/3 



9/0 

0.28 

42.29 


Number of Planck detections. (Z,) Number of Planck upper limits. 

* The number of detections is too low to allow a reliable estimation of p. 
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Fig. 12. cthf vs. u\ y diagram for v Break = 70 GHz for different 
blazar classes: Flat Spectrum Radio Quasars (FSRQ), BL Lac 
objects (BL Lac), and AGN of uncertain classification. 


the time of the Planck observations. From Table |23] we see that, 
although the significance of the correlation in the various sam- 
ples and blazar classes substantially increases, it is never very 
high and the scatter remains. 

This result may be due to the limited flux dynamic range in 
the y-ray band in the Fermi-LAT sample and to the small number 
of objects (for the case of BL Lacs). However, the low level of 
correlation could occur if the microwave flux density represents 
the superposition of multiple synchrotron components while the 
simultaneous y-ray flux represents the emission from a single 
dominant component which may be active only for a short time. 
For this reason in the bottom panel of Fig.[l4j we plot the Planck 
ERCSC flux densities (which in most cases represents the flux 
density averaged over more than one Planck observation) versus 
the Fermi - LAT flux averaged over the entire 27 month period. 
In this case the correlation is highly significant, however, the 
many upper limits in both energy bands clearly imply that the 
dispersion is very large. 

7.3. X-ray vs. y-ray 

Figure |T5] plots the 2 month y-ray flux versus the simultaneous 
X-ray flux for all sources observed by Swift that do not show sig- 
natures of thermal contamination in the X-ray spectrum. Open 



Fig. 13. The Planck 143 GHz flux density plotted vs. the simul- 
taneous Swift XRT X-ray flux for all sources where the X-ray 
flux is expected to be due to the inverse-Compton component 
are plotted. Blazars with significant X-ray contamination due to 
accretion have been excluded from the plot. 


red circles represent HSP sources, where the X-ray flux is due 
to the tail of the synchrotron emission, while black filled circles 
are LSP and ISP sources for which the X-ray flux is related to 
the inverse-Compton radiation. This distinction is needed in or- 
der to properly compare sources where the emission is due to the 
same mechanism. We therefore compute the y-ray vs. X-ray cor- 
relation coefficient only for LSP and ISP sources. The Spearman 
rank test shows moderate evidence for a correlation which is less 
significant in the longer integrations: P = 2.55% for simultane- 
ous data, P = 6.9% for 2 months integration, or P = 17.4% for 
27 months integration. 

Figure[l6]plots the power law spectral index in the Swift XRT 
energy band (0.3-10 keV) versus the spectral slope in the Fermi- 
LAT y-ray band (lOOMeV-lOOGeV) derived using the 2-month 
y-ray data. A significant correlation is present, confirming the 


shows the same plot built with simultaneous y-ray data; although 
the number of objects is lower the correlation is still present. 


results of Abdo et al. ([2010a I. The top right corner of Fig. 16 


8. Discussion and conclusions 

We have collected simultaneous Planck, Swift, Fermi-LAI, and 
ground-based multi-frequency data for 105 blazars included in 
three statistically well defined samples characterized by flux lim- 
its in the soft X-ray (0.1-2.4keV, ROSAT ), hard X-ray (15- 
150 keV, Swift- BAT), and y-ray (E > 100 Me V, Fermi- LAT) en- 
ergy bands. This study complements a similar study of a sample 
of 104 radio-bright AGN (JyjGHz > lJy) (Planck Collaboration 
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Table 23. Microwave (143 GHz) vs. y-ray flux correlation parameters. 


y-ray 

integration 

period 

Class 

Fermi- LAT sample 
No. of p P(%) 

sources 
(a)/(b)/(c) 

Swift- BAT sample 
No. of p P(%) 

sources 
(a)/(b)/(c) 

RO S AT /R ASS sample 
No. of p P(%) 

sources 
(a)/(b)/(c) 

Radio sample 
No. of p 

sources 
(a)/(b)/(c) 

P(%) 


All 

16/5/32 

0.22 

12.4 

8/6/20 

0.42 

2.49 

6/20/25 

0.42 

1.42 

22/2/71 

0.35 

0.07 

Simult. 

FSRQs 

10/3/16 

0.23 

23.4 

5/2/12 

0.39 

9.92 

4/3/6 

0.65 

3.97 

15/1/46 

0.33 

1.19 


BL Lacs 

5/2/9 

0.55 

4.76 

2/2/2 

...* 


2/11/4 



4/0/10 

0.38 

17.6 


All 

36/5/9 

0.29 

4.18 

11/7/20 

0.61 

0.07 

8/21/18 

0.31 

6.00 

37/3/56 

0.40 

0.01 

2 months 

FSRQs 

19/3/6 

0.17 

38.7 

6/3/15 

0.59 

0.89 

4/4/9 

0.51 

6.83 

23/2/42 

0.36 

0.39 


BL Lacs 

11/2/1 

0.72 

0.92 

3/3/2 



3/11/5 



10/0/6 

0.42 

10.2 


All 

47/3/0 

0.47 

0.10 

22/16/7 

0.67 

<0.01 

12/23/18 

0.40 

0.99 

75/4/16 

0.48 

<0.01 

27 months 

FSRQs 

25/3/0 

0.48 

1.27 

16/3/4 

0.54 

1.59 

8/3/7 

0.85 

0.15 

47/2/19 

0.48 

0.01 


BL Lacs 

14/0/0 

0.72 

0.91 

3/4/0 



2/14/2 



16/0/0 

0.55 

3.22 


(a) Number of sources detected both by Planck and Fermi-LAT. (i) Number of Planck upper limits. <c) Number of Fermi- LAT upper limits. 
* The number of detections is too low to allow a reliable estimation of p. 
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Fig. 14. Top panel: The Planck flux density at 143 GHz is plot- 
ted against the simultaneous Fermi-LAI flux. Middle panel : The 
Planck flux density at 143 GHz is plotted against the Fermi-LAI 
flux integrated over the 2 month period entered on the Planck- 
Swift observations. Bottom panel: The Planck ERCSC flux den- 
sity at 143 GHz (flux averaged over more than one Planck sur- 
vey) is plotted versus the y-ray fluxes averaged over the entire 
27 month Fermi-LAI observing period. 


201 le) . Altogether the four samples include a total of 175 
distinct objects. The acquisition of this unprecedented multi- 
frequency/multi-satellite data set was possible thanks to coop- 



Fermi 2 Month Flux (ph cm-V) 


Fig. 15. The Swift XRT flux is plotted against the 2 month Fermi- 
LAI flux for the sources included in the three high energy flux 
limited samples. Open red circles represent HSP sources, i.e., 
high synchrotron peaked BL Lacs where the X-ray flux is dom- 
inated by synchrotron radiation. Filled black circles represent 
LSP sources, i.e., blazars with low synchrotron peak, for which 
the X-ray flux is dominated by inverse-Compton radiation. 


eration between the Planck, Swift and Fermi- LAT teams, who 
agreed to share data and to organize an extensive program 
of multi -frequency observations involving over 160 Swift ToO 
pointings scheduled when the blazars were within the field of 
view of the Planck instruments. 

We have used this unique multi-frequency dataset to build 
well-sampled, simultaneous SEDs of all the blazars included in 
our high-energy selected samples. This collection of SEDs is 


an improvement over previous compilations (e.g., Abdo et al. 
2010a) because: a) the SEDs presented here are strictly simul 


taneous while in |Abdo et~ak1 ( |20 lOaj) the multi-frequency data 
were collected over a period of up to 9 months centered on the 
first three months of Fermi-LAI operations; b) the sources were 
selected according to different statistical criteria allowing us to 
probe the blazar parameter space from widely different view- 
points; and c) we took care to identify and separate radiation 
components unrelated to the emission from the jet, such as the 
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Photon index (y-ray band) 


Fig. 16. The power law spectral index in the Swift XRT energy 
band (0.3-10 keV) is plotted versus the slope in the Fermi-LAT 
y-ray band (lOOMeV-lOOGeV) derived using Fermi-LAT data 
integrated over a period of 2 months. The same plot built with si- 
multaneous y-ray data is shown in the inset at the upper right. A 
clear anti-correlation is present; the Spearman test gives a prob- 
ability of less than 0.01% that the correlation is due to chance, 
even in the simultaneous data. 


light from the blazar’s host galaxy and the radiation due to the 
accretion onto the central black hole, which often contaminate 
the non-thermal blazar spectrum in the optical, UV, and X-ray 
bands. 


Our findings are broadly consistent with those of Abdo et al. 
( 201 0a| . However, our larger samples selected in different parts 
of the electromagnetic spectrum, better wavelength coverage, 
level of simultaneity, and separation of the emission compo- 
nents, allowed us to make significant progress in several areas. 
The use of four widely different samples has allowed us to in- 
vestigate the consequences of selection effects on the estimation 
of important parameters such as (v* eak >, average Compton dom- 
inance, etc., especially for the case of BL Lacs. 

Some of our sources have been observed simultaneously by 
Swift, Planck, and Fermi-LAT during more than one Planck sur- 
vey. In these cases here we have presented only the data collected 
during the first observation. Multiple simultaneous observations 
of a subset of our blazars and detailed model fitting of the SEDs 
will be the subject of future papers. The main results of this work 
are discussed below. 


8.1. Fermi-LAT detection statistics and the effects of 
variability 

The percentage of Fermi-LAT detected sources during the simul- 
taneous integrations, typically lasting about one week, ranges 
from < 40% in the Fermi-LAT sample to 20-25% in the radio 
and X-ray selected samples (see Table m When 2-month in- 
tegrations centered on the Planck-Swift observations are consid- 
ered, these percentages grow to 80% in the Fermi-LKL sample 
and to ~ 35% in the other samples. But even with the 27-month 


Fermi-LALL integrations available at the time of writing, many of 
the blazars belonging to the radio, soft and hard X-ray selected 
samples remain undetected. 

It is interesting to note that the detection rate is quite differ- 
ent in FSRQs and BL Lacs: if we exclude the Fermi-LKL sample 
where all the objects have been detected (by definition), the per- 
centage of detections is 95% for BL Lacs and only ~60% for 
FSRQs, with values ranging from 72% in the radio sample to 
just 53% in the RASS sample (see Table [12}. 

A comparison of our simultaneous data with published and 
archival measurements shows that the use of non-simultaneous 
data in the SED of blazars typically introduces a scatter of about 
a factor 2 in the microwave band, and of a factor of up to ten or 
more in the X-ray and y-ray bands. 


8.2. The spectral slope in the radio-sub-millimetre region 


We confirm that the energy spectrum of blazars in the radio- 
microwave spectral region is quite flat, with an average slope 
(a) ~ 0 (/(v) oc v a ) up to about 70 GHz, above which it steep- 
ens to (a) ~ -0.6. This behaviour is very similar to that ob- 
served in the sample of radio bright blazars considered by |Planck 
[Colla boration (201 le ! (see also |Tucci et al.|20TT) confirming the 
findings of|Abdo et al.|(|2010a[) that the radio to microwave part 


of the spectrum is approximately the same in all blazars (FSRQs 
and BL Lacs) independently of the selection band. However, the 
spectral slope of BL Lacs above ~70 GHz is flatter than that of 
FSRQs with (q-hf) = -0.51 + 0.07 for BL Lacs compared to 
<a HF > = -0.73 ± 0.04 for FSRQs (see Table [TF] and Fig. [Tl}. 
A KS test, performed on the subsamples of FSRQs and BL Lacs 
with radio flux density larger than 1 Jy, gives a probability of less 
than 3% that the two samples come from the same parent popu- 
lation. This difference in the high frequency spectral index may 
be due to a primary dependence of othf on Vp eak , as illustrated in 

Fig. jlTJ because only BL Lacs show high Vp eak values. 


8.3. Synchrotron self-absorption 

We searched for signatures of synchrotron self-absorption, 
which in simple homogeneous SSC models is expected to cause 
strong spectral flattening below ~ 100 GHz, but we found no ev- 
idence for common behaviour; indeed, the average spectrum in 
that region steepens instead of flattening. Possible cases where 
some evidence of synchrotron self-absorption may be present are 
PKS 0454-234, PKS 0521-36 (Fig. [28}, S4 0917+44 (Fig. [30}. 
PKS 1127-145 (Fig. [32} and 3C 454.3 (Fig. [41} . 

8.4. Non-thermal versus disk radiation 

In several blazars the optical/UV light is significantly contam- 
inated by thermal/disk radiation (known as the blue bump, see 
Figs. [6] and [7} while the soft X-ray flux is contaminated by ra- 
diation produced in the accretion process in approximately 25% 
of the blazars in our samples (see Fig. [9}. In some of the closest 
sources the optical light is instead contaminated or dominated 
by the emission from the host galaxy. Ignoring this contamina- 
tion may cause an overestimate of the position of both v s . and 

J r peak 

Vp aak by 0.5 dex or more. 

We investigated the relationship between the radiation due 
to accretion and to the jet, which co-exist in most FSRQs, us- 
ing the parameter ctr onhermain defined as the spectral slope be- 
tween the 5 GHz radio (non-thermal) flux and the 5000 A op- 
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Fig. 17. The spectral slope cthf is plotted against v peak for all 
the sources in our samples with a measured value of q-hf and 
S5 ghz > 1 Jy. High v s eak blazars tend to be harder microwave 
sources. 
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Fig. 18. The rest-frame synchrotron v peak distributions of FSRQs 
and BL Lacs in different samples. Black solid histograms: radio 
sample; red dashed histograms: Fermi - LAT sample; green dot- 
dashed histograms: Swift BAT sample; blue dotted histograms: 
RASS sample. 


tical flux that can be attributed to the blue-bump/disk emis- 
sion. In the blazar paradigm this quantity depends both on the 
amount of relativistic amplification (which only affects the non- 
thermal radiation from the jet) and the intrinsic ratio between 
non-thermal/jet radiation and disk emission. We estimated the 
value of aR-o(Thermai) in all sources for which we have optical 
data of good quality. Since BL Lacs do not show broad lines 
only lower limits to crR-ocrhermai) can be derived. 

Figure [8] shows the distribution of «'r o(Thermai) for all our 
samples, which range from ~0.4 to just above 1.0 with peak 
values between 0.6 and 0.8. The distributions are all similar, 
with the largest difference being between the Fermi-L AT and 
the RASS samples (a KS test gives a probability of 3.6% that the 
two distributions are from the same parent population), possibly 
reflecting differences in the amplification factor or a different ra- 
tio of accretion to jet emission in y-ray and soft X-ray selected 
blazars. We note that the distribution of the arR_o(Thermai) limits 
for BL Lacs is consistent with intrinsic values of trR-o(Thermai) for 
BL Lacs that are within the range of values observed in FSRQs. 


8.5. The distribution of rest-frame peak energies 


The distribution of rest-frame synchrotron peak energies ( v 'p Cak 1 
of FSRQs is very similar in all our samples with an average of 
( v peak^ ~ 10 131±al Hz and a dispersion of only 0.5 dex. But for 
BL Lacs the value of (Vp eak ) is at least one order of magnitude 
larger than that of FSRQs, the exact value depending consider- 
ably on the selection method (see Table[20|and Fig.[T8j). The dis- 

■ivak 


tributions of vj£,, for FSRQs and BL Lacs also differ, but not as 


much as those of v 3 e . ik (see Fig. 19}. The majority of the sources 


in all the samples (both FSRQs and BL Lacs) peak between 
10 21 and 10 23 Hz, with a few extreme HSP BL Lacs reaching 
~ 10 26 Hz. 
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Fig. 19. The rest-frame inverse-Compton v pL , ak distributions of 
FSRQs and BL Lacs in different samples. Black solid his- 
tograms: radio sample; red dashed histograms: Fermi - LAT sam- 
ple; green dot-dashed histograms: Swift BAT sample; blue dotted 
histograms: RASS sample. 


8.6. Correlations between fluxes and other blazar 
parameters 

Despite the strict simultaneity of our data, plots of fluxes in dif- 
ferent spectral regions (microwave vs. X-ray, microwave vs. y- 
ray. X-ray vs. y-ray) still show a large scatter (see Figs. 13 14 


tween the radio and y - ray have been reported (e.g. Kovalev et al. 

2009 

Leon-Tavares et al. 2010; 

Giroletti et al. 2010; Ghirlanda 

et al. 

2010; Mahony et al. 2010 

Peel et al.||2010; Linford et al. 

2011 

Ackermann et al.]201 1). T 

te difference might be due to the 
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different synchrotron peak energies of the objects in the samples. 
Even in simple SSC scenarios this introduce a scatter in the cor- 
relation between the fluxes (e.g., for the same radio flux, an ob- 
ject with larger Vp eak is expected to produce more y-rays than one 

with a smaller Vp eak ). The large scatter present in Fig. |l4| could 
also imply that y-ray emission is due to components not always 
directly related to radiation in other energy bands, e.g., multiple 
SSC components (see also |Abdo et "alTjl201 0a| . A good corre- 
lation is present between the X-ray and y-ray spectral slopes. 


however (see Fig. 16 1 . 


We confirm the correlations between the Fermi-L AT spec- 
tral index and the SED peak energies v s , and v ,c , found by 

r ° peak peak J 

[Abdo et ak] ( [2010a| . As an illustration of the agreement, Fig.|20| 
plots the Fermi - LAT spectral slope estimated using the full 27- 
month data set as a function of vj^ ak . The gray points represent 
the y-ray spectral slopes estimated using the quasi simultane- 
ous 2-month dataset. These points, plotted without the much 
larger error bars to avoid confusion, clearly cluster around the 
27-month data, confirming the correlation. The solid line repre- 
sents the best fit obtained by Abdo et al. (2010a|>. 


where y** ak is the Lorentz factor of the electrons radiating at the 
peak energy. This is related to the observed peak frequency of 
the observed photon spectrum by: 


Tpeak K 


peak 

~B5 


1/2 


(4) 


where v* eak is the synchrotron peak frequency in the rest-frame 
of the emitting region, B is the magnetic field, and 5 is the usual 
Doppler factor (e.g. Urry & Padovani |1995[ >. 

In objects where v^ c; k is larger than ~ 10 15 Hz the Thomson 
approximation is no longer valid and the inverse-Compton scat- 
tering occurs under the Klein-Nishina (KN) regime. By means 
of Monte Carlo simulations, |Abdo et ah] ( |2010a] ) estimated the 
area covered by SSC models in the plane l°g(t / p eak )-l°g(Tp e S ak )- 
This area is delimited by the solid contour lines shown in Fig. 21 
where we plot the log(y s ^) of our sources, calculated from the 

|3j versus the rest- 

peak 


observed values of Vp eak and vj^’ ak using eq. 
frame log(vL t ). 



Fig. 20. The y-ray spectral index, estimated from the entire 27- 
month Fermi - LAT data set of all sources in our samples is plotted 
against log(y I ^ ak ). Gray points, plotted without statistical errors 
to avoid confusion, represent the Fermi -LAI spectral slopes es- 
timated using the 2-month integrations. The solid line is the best 
fit given by Abdo et al. ( 20 1 Qa[> 
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Fig. 21. log(yS^), obtained from Eq. |3j for all the objects in 
our samples for which we could obtain Vp eak and vj^ ak is plot- 
ted against log(V^ ak ) in the rest frame of the blazars. The two 
black lines delimit the area predicted by simple homogeneous 
SSC models obtained through extensive Monte Carlo simula- 
tions (see|Abdo et al.|2010a| for details). 


8. 7. Comparison with the expectation of simple SSC models 

As discussed in Abdo et al. ( |2010a| >, simple SSC models pre- 
dict that in the Thomson regime the peak frequency of the syn- 
chrotron (Vp eak ) and inverse-Compton (v' aak ) components are re- 
lated by: 


,ic 

peak 


peak 



(3) 


As in the case of the bright y-ray blazars considered by |Abdo| 
et al. (2010a|>, only a few objects are inside or close to the SSC 


area, implying that simple SSC models cannot explain the SED 
of many blazars in our samples. This conclusion is supported 
by other evidence including the lack of a strong correlation be- 
tween radio and y-ray and the fact that the local optical/UV 
spectral slope in some objects, although not contaminated by 
thermal radiation, is inconsistent with the expectations of SSC. 


See, e.g., PKS 0426-380 (J0428-3756, Fig. [28}, PKS 0851+202 
(J0854+2006, Fig. [30}, S5 1803+784 (J1800+7828, Fig.[38}. 
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Fig. 22. The logarithm of the Compton dominance is plotted as 
a function of log(v s , ) for all sources for which v s , and v IC , 

ov peak 7 peak peak 

could be reliably determined. 


Note, however, that the 0 f blazars that were not detected 

by Fermi-LAI and for which we could obtain a limit to v IC . 

J peak 

(30 to 40% of the FSRQs in the radio and X-ray selected sam- 
ples; see Fig. 10 for one example) are plotted as upper limits in 
Fig. [2TJ many of these limits are close to or inside the SSC area 
and therefore the SEDs of these objects are likely consistent with 
simple SSC emission. 


8.8. The Compton dominance of blazars 

The Compton dominance (CD, defined as the ratio between the 
inverse-Compton and synchrotron peak luminosities) is a crucial 
parameter for the study of blazar physics, as it is strictly related 
to the location of the maximum power output in the energy spec- 
trum of a blazar. 

Figure [22]plots the CD values, estimated from our SEDs, as 
a function of Vp eak , showing that log(CD) ranges from about -0.5 
to about 2. The larger values are always associated with LSP ob- 
jects while HSP sources always show log(CD) lower than ~ 0.5. 
In this figure the two blazar subclasses appear to be quite dif- 
ferent, with BL Lacs having significantly lower CD values, even 
when their v s . values are equal to those of FSRQs. To better un- 

peak 1 ^ 

derstand this difference, in Fig.[23]we plot the CD distributions 
of FSROs and BL Lacs for different samples and v s , intervals. 

^ r peak 

FSRQs included in the Fermi-LAI sample, which are y- 
ray bright by definition, show a CD distribution peaking at large 
values. Note that in this sample we also applied a radio flux- 
density limit of 1 Jy, so the sources below the radio cut must be 
on average more Compton dominated than those in our sample. 
This implies that the distribution of purely y-ray selected blazars 
must be even more strongly peaked at high CD values than that 
of the Fermi-LAT sample. Considering instead FSRQs selected 
in the radio and in the X-ray bands, we get quite a different pic- 
ture, with a broader distribution extending to values of less than 
1. Moreover, about 30 — 45% of FSRQs in the radio, soft X-ray, 
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Fig. 23. Distributions of the Compton dominance for FSRQs and 
BL Lacs of the LSP and HSP type. A significant fraction of 
FSRQs in the radio, RASS, and BAT samples have not been de- 
tected in the y-ray band and therefore only limits to the Compton 
dominance (shown as a dashed histogram; only the radio sample 
to avoid confusion) can be calculated. The large difference be- 
tween the CD distribution of FSRQs in the Fermi-LAT and in 
the other samples illustrates the strong bias that y-ray selection 
induce. 


and hard X-ray selected samples are not detected by Fermi-LAT 
and therefore they must populate the part of the CD distribution 
with low CD values. This is shown by the dotted red histogram, 
which also includes upper limits to the CD (estimated as the ratio 
of the upper limit to v ,c . F(v ic . ) and v s , F(v s , ) where limits 

l ^ peak v peak 7 peak x peak 7 

to v ,l( . and v IC . F(v ic , ) are obtained by fitting the X-ray data 

peak peak v peak 7 jo j 

together with the 27-month Fermi-LAI upper limits as shown in 
Fig. [10] 


8.9. The blazar sequence 


The top panel of Fig. [24] plots the logarithm of the bolometric 
power, represented by the sum of the synchrotron and inverse- 
Compton peak luminosities [L Bo i ~ v® eak T(v^ eak ) + vj^Lfv^)] 


as a function of log(v peak ) for all sources in the four samples 


considered in this paper and in (IPlanck Collaboration||201 leb 

„ ... . „ e ...... ... 


for which an estimate of v s , and of the bolometric luminos- 

peak 

ity was possible. We use this plot to test the relationship known 
as the Blazar Sequence, that is the strong anti-correlation be- 
tween bolometric luminosity and v s , claimed by Fossati et al. 
(1998) and Ghisellini et al. (1998) which has been the subject of 


lively debate (e.g., Giommi et al. |1999| Padovani et al. 2003; 


Caccianiga & Marcha 2004| Nieppola et al.||2006| |Padovani| 
2007[|Ghisellini & Tavecchio|2008[ ). 


We stress that in order to properly test for the existence of 
such a relationship it is mandatory to use samples that are unbi- 
ased, that is selected in such a way so as not to favour any part 
of the CBoi-fpe^ diagram. Although our samples are statistically 
well defined they are not unbiased from this viewpoint for the 
following reasons: 


22 


P. Giommi et al. : Simultaneous Planck, Swift, and Fermi observations of blazars 



12 ' 14 ' 16 ' W 


Log(v s peak , Hz) 


Fig. 24. Top panel: The bolometric luminosity (represented by 
the sum of synchrotron and inverse-Compton peak luminosities) 
is plotted against y s eak . The non-thermal optical light of blazars 
with no strong emission lines above the dashed line would be 
bright enough to swamp the emission from the host galaxy 
making the source appear featureless thus hampering any red- 
shift measurement. The upper limits represent BL Lacs with no 
known redshift are estimated assuming that the non-thermal light 
is larger than ten times that of a giant elliptical host galaxy (see 
text for details). The bottom panel shows the same plot for the 


subsample of sources satisfying the selection criteria of Fossati 
|etak| ( fT998| l. 


redshift measurement ( |Massaro et al.||2010[ |Shaw et aki r2009 


2010). All the BL Lacs in our samples for which redshifts are 


unknown are plotted in Fig. [24]as lower limits estimated assum- 
ing that their non-thermal light is larger than 10 times the op- 
tical light of the host galaxy. The fact that some BL Lacs with 
known luminosity (red points in Fig. [24| are above the red line 
is because their redshift was measured from emission lines with 
rest-frame equivalent width below the 5A limit; they might be 
intermediate objects between BL Lacs and FSRQs (see. e.g., 
[Ghisellini et al.[201 1 [ ). 

c) The different radio flux density cuts applied to our sam- 
ples imply that different subsamples probe different parts of the 
radio luminosity function. The radio and y- ray selected samples, 
which are defined using a high radio flux limit (S > 1 Jy), probe 
the high-luminosity end of the radio luminosity function, which, 
for Vp eak — 10 13 Hz and the observed redshift and CD distribu- 
tions, translates into L Bol > 10 46 erg s 1 . The soft and hard X-ray 
selected samples have a radio flux density cut of S - 0. 1-0.2 Jy, 
or about one order of magnitude fainter than that of the radio and 
Fermi - LAT samples, and therefore can include significantly less 
powerful sources, down to bolometric luminosities of the order 
of 10 44 ergs~ ! . 

If we remove the v* eak dependency on the selection method 
by considering FSRQs only (which show the same v s eak distri- 
bution in all samples), we see that the luminosity values span 


ergs 

10 125 


about five orders of magnitude from ~ 10 44 to ~ 1 0 49 
and show no trend with Vp eak , which ranges between ~ 

and ~ 10 14 Hz. No obvious correlation is present in each sample 
separately or in the union of the four samples. 

Finally we consider the case of the subsample of sources that 
satisfy the same conditions as in Fossati et al. (1998), that is 
S 5 ghz > 2 Jy for FSRQs, S 5 ghz > 1 Jy for radio selected BL 
Lacs, no restrictions for X-ray selected BL Lacs, and excluding 
all BL Lacs with no redshift information. This case is illustrated 
in Fig. [24] (bottom panel) which shows a trend very similar to 
that presented in Fossati et al. ( 1998). 

Taking take into account all of the above, we conclude that 
our data do not show a correlation of the type predicted by 
Fossati et al. ( 1998[ >, owing to the presence of low luminosity 


LSP objects and the difficulty of measuring redshifts for likely 
high-luminosity HSP sources. That such a correlation becomes 


evident when the objects are selected by the criteria of Fossati 
et al. (|1998|l supports the hypothesis that the correlation might 


be the result of a selection effect. 


a) As shown in Figs. 


18 


and 


19 


the distribution of 


peak 


strongly depends on the selection method; hence, depending on 
the mix of the samples, we obtain a different sampling of the 
parameter v® eak . 


b) The area above the dashed line, which represents the lu- 
minosity above which the non-thermal emission of a blazar com- 
pletely dominates the observed optical fluxj^] cannot be pop- 
ulated by blazars with no emission lines, like BL Lacs, as in 
this case they would appear completely featureless and therefore 
their redshift could not be measured. In this respect we note that 
over 40% of the BL Lacs in the BZCAT catalog, and an even 
larger fraction of the Fermi-LAI detected BL Lacs, still lack any 


12 We assume a non-thermal luminosity one order of magnitude larger 
than that of the host giant elliptical galaxy of luminosity equal to that 
found by |Scarpa et al.|p000| > and by Urry et al7| i }2000[ > and verified by 
us to fit our SEDs. 


8. 1 0. Selection effects and sample composition 


Our choice of considering flux-limited samples selected in 
widely different parts of the electromagnetic spectrum (radio, 
soft X-ray, hard X-ray, and y-ray) allowed us to emphasize the 
strong selection biases that can affect important physical pa- 
rameters such as the peak energy of both the synchrotron and 
inverse-Compton components (see Figs. 18 and [19} and the 
Compton dominance (see Fig. |23) . Since FSRQs and BL Lacs 
show significantly different v 3 cak distributions, these selection bi- 
ases also strongly affect the composition of the samples in terms 
of relative abundance of blazar subclasses (FSRQs vs. BL Lacs, 
LSPs vs. HSPs), as is apparent from Table[4] 

Radio-selected samples include sources that are bright 
in the radio band. If there is no correlation between radio 
flux/luminosity and other parameters such as v 3 aak , Compton 
dominance, etc., this is the best selection for measuring the dis- 
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tributions of these important physical parameters. If instead there 
is a strong correlation between radio luminosity and v® eak , then 

the distribution of Vp eak should strongly depend on the radio flux 
limit. 

X-ray selection favors high v* eak (and consequently high 
v peak) sources which are much brighter at X-ray frequencies than 
low v s . for the same radio flux. X-ray flux limited samples 

peak J r 

are therefore much richer in high Vp eak BL Lacs (HBLs or HSP 
sources) than radio-selected samples. This selection effect has 
been well known ever since the first soft X-ray surveys became 
available. 

y-ray flux limits favor bright y- ray objects and there- 
fore highly Compton-dominated sources. Fermi-L AT TS-limited 
samples favor sources with flat y-ray spectral slopes or high 
v IC , sources. This explains the overabundance of HSP blazars 

peak r 

(only BL Lacs) and of high CD blazars (only FSRQs) in Fermi- 
LAT catalogs. 
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J0010+1058 - III ZW 2 


J0017+8135 - S5 0014+813 








Fig. 25. The SED of IIIZW2 (J0010+1058, top left), S5 0014+813 (J0017+8135, top right), 1ES 0033+595 (J0035+5950, middle 
left), Mkn 348 (J0048+3157, middle right), lJy 01 18-272 (J0120-2701, bottom left), and S40133+47 (J01 36+4751, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J0204-1 701 - PKS 0202-17 


J021 0—51 01 - PKS 0208-512 








Fig. 26. The SED of PKS0202-17 (J0204-1701, top left), PKS0208-512 (J0210-5101, top right), GB6J 0214+5145 
(J0214+5144, middle left), PKS 0215+015 (J0217+0144, middle right), lJy 0212+735 (J0217+7349, bottom left), and 
lJy 0218+357 (J0221+3556, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data inte- 
grated over 2 months, Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated 
over 27 months are shown in blue; literature or archival data are shown in light gray. 
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Fig. 27. The SED of 4C 28.07 (J0237+2848, top left), PKS 0235+164 (J0238+1636, top right), NGC 1275 (J0319+4130, middle 
left), PKS 0332-403 (J0334-4008, middle right), NRAO 140 (J0336+3218, bottom left) and PKS 0420-01 (J0423-0120, bottom 
right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and 
non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature 
or archival data are shown in light gray. 
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J0428-3756 - PKS 0426-380 


J0433+0521 - 3C 120 








Fig. 28. The SED of PKS 0426-380 (J0428-3756, top left), 3C 120 (J0433+0521, top right), PKS 0454-234 (J0457-2324, middle 
left), PKS 0521-36 (J0522-3627, middle right), PKS 0528+134 (J0530+1331, bottom left), and PKS 0537-441 (J0538-4405, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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J0539-2839 - Uy 0537-286 


J0550-3216 - PKS0548-322 








Fig. 29. The SED of lJy 0537-286 (J0539-2839. top left), PKS 0548-322 (J0550-3216, top right), IR AS -L 06229-643 
(J0623-6436, middle left), PKS 0735+17 (J0738+1742, middle right), B2.2 0743+25 (J0746+2549, bottom left), and S4 0814+425 
(J0818+4222, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, 
Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are 
shown in blue; literature or archival data are shown in light gray. 
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J0824+5552 - OJ 535 


J0841+7053 - 4C71.07 








Fig. 30. The SED of OJ535 (J0824+5552, top left), 4C71.07 (J0841+7053, top right), PKS 0851+202 (J0854+2006, middle left), 
B2 0912+29 (J0915+2933, middle right), S4 0917+44 (J0920+4441, bottom left), and PKS 0921-213 (J0923-2135, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J0957+5522 - 4C 55.17 


J1 01 5+4926 - 1 H 1013+498 








Fig. 31. The SED of 4C55.17 (J0957+5522, top left), 1H 1013+498 (J1015+4926, top right), 1RXSJ 105837.5+562816 
(J1058+5628, middle left), 4C01.28 (J1058+0133, middle right), PKS 1057-79 (J1058-8003, bottom left), and Mkn421 
(J 1 104+3812, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, 
Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are 
shown in blue; literature or archival data are shown in light gray. 
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J1 127-1 857 - PKS 1124-186 


J1 130-1449 - PKS 1127-145 








Fig. 32. The SED of PKS 1 124-186 (J1127-1857, top left), PKS 1127-145 (J1 130-1449, top right), B2 1128+31 (J1131+3114, 
middle left), S5 1133+704 (J1136+7009, middle right), PKS 1144-379 (J1 147-3812, bottom left), and 4C 49.22 (J1 153+4931, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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J1 159+2914 - 4C 29.45 
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Fig. 33. The SED of 4C29.45 (J1 159+2914. top left), ON 325 (J1217+3007, top right), PKS 1217+02 (J1220+0203, middle 
left), ON 231 (J 1221+2813, middle right), PKS 1219+04 (J1222+0413, bottom left), and 3C273 (J1229+0203, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J1 246-2547 - PKS 1244-255 


J1 248+5820 - PG 1246+586 








Fig. 34. The SED of PKS 1244-255 (J1246-2547, top left), PG 1246+586 (J1248+5820, top right). 3C279 (J1256-0547, middle 
left), lJy 1302-102 (J1305-1033, middle right), lJy 1308+326 (J1310+3220, bottom left), and GB6B 1347+0955 (J1350+0940, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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J1 407-2701 - 1 WGAJ 1407.5-2700 


J1 41 9+0628 - 3C 298.0 








Fig. 35. The SED of 1WGAJ 1407.5-2700 (J 1407-2701, top left) and of 3C298.0 (J1419+0628, top right), CS0643 
(J1423+5055, middle left), of PG 1424+240 (J1427+2348, middle right), 1RXSJ 145603.4+504825 (J1456+5048, bottom left), 
and PKS 1502+106 (J1504+1029, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data inte- 
grated over 2 months, Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated 
over 27 months are shown in blue; literature or archival data are shown in light gray. 
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J1 507+041 5 - BZQJ 1507+0415 


J1 51 0-0543 - 4C -05.64 








Fig. 36. The SED of BZQJ 1507+0415 CJ 1507+0415. top left), 4C -05.64 (J1510-0543, top right), APLib (J1517-2422, middle 
left), PG 1553+113 (J1555+1 1 11, middle right), WE 1601 + 16W3 (J1603+1553, bottom left), and OS -237.8 (J1625-2527, bottom 
right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and 
non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature 
or archival data are shown in light gray. 
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J 1635+3808 - 4C 38.41 


J 1640+3946 - NRA0 512 








Fig. 37. The SED of 4C38.41 (J1635+3808, top left), NRA0512 (J1640+3946, top right), 3C345 (J1642+3948, middle left), 
Mkn501 (J1653+3945, middle right), ARP 102B (J1719+4858, bottom left), and 1ES 1741 + 196 (J1743+1935, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J1 751 +0939 - OT 081 


J1 800+7828 - S5 1803+784 








Fig. 38. The SED of OT081 (J1751+0939, top left), S5 1803+784 (J1800+782, top right), PKSB 1830-210 (J1833-2103, middle 
left), PKS 1833-77 (J1840-7709, middle right), 2E 1908.2-201 (J 19 11 -2006, bottom left), and PMNJ 1923-2104 (J1923-2104, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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J1 924-291 4 - OV -236 


J 1959+6508 - 1ES 1959+650 








Fig. 39. The SED of OV-236 (J1924-2914, top left), 1ES 1959+650 (J1959+6508, top right), lJy 2005-489 (J2009-4849, middle 
left), PKS 2052-47 (J2056-4714, middle right), lJy 2126-158 (J2129-1538, bottom left), and S3 2141 + 17 (J2143+1743, bottom 
right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and 
non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature 
or archival data are shown in light gray. 
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J21 47+0929 - Uy 2144+092 


J21 48+0657 - 4C 06.69 








Fig. 40. The SED of lJy2144+092 (J2147+0929, top left), 4C 06.69 (J2148+0657, top right), PKS2149-307 (J2151-3027, mid- 
dle left), BLLac (J2202+4216, middle right), 4C 31.63 (J2203+3145, botton left), and PKS 2204-54 (J2207-5346, botton right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J2209-4710 - NGC7213 


J2229-0832 - PKS 2227-08 
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Fig. 41. The SED of NGC 7213 (J2209-4710, top left), PKS 2227-08 (J2229-0832, top right), PKS 2227-399 (J2230-3942, mid- 
dle left), 4C 11.69 (J2232+ 1143, middle right), 3C 454.3 (J2253+ 1608, bottom left), and PKS 2300-18 (J2303-1841, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J2327+0940 - PKS 2325+093 


J2333-2343 - PKS 2331 -240 




J2347+5142 - 1 ES 2344+51 4 



Fig. 42. The SED of PKS 2325+093 (J2327+0940. top left), PKS 2331-240 (J2333-2343, top right), and 1ES 2344+514 
(J2347+5142, middle). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, 
Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are 
shown in blue; literature or archival data are shown in light gray. 
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Table 6. Planck data integrated over a single survey during the Swift observations. 
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Columns 1 and 2: source name; Columns 3 and 4: observation start and end times; Columns 5-13: flux densities in units of Jy at 30, 44, 70, 100, 143, 217, 353, 545, and 857 GHz. 



Table 7. Swift - UVOT data. 
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Table 7. Swi/i-UVOT data. Continued. 
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1058-8003 PKS 1057-79 2010-08-30 -2.00 0.63 1.0 0.9 0.08 

1635+3808 4C 38.41 2010-03-07 -2.00 Oil L3 U 0.25 

Column 1 and 2: source name; Column 3: Swift observation date; Column 4: best fit photon index r t ; Column 5: Galactic N^; 
Column 6 and 7: 0. 1-2.4 and 2-10keV X-ray fluxes; Column 8: reduced^ 2 ; Column 9: number of degrees of freedom. 
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Table 13. Fermi-LK T data at 2 bins/decade simultaneous with the Planck observations. 
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Table 14. Fermi-LAT data at 1 bin/decade simultaneous with the Planck observations. 


J2000.0 name 

Source name 

Flux density (10 u 'phcm 2 s 

MeV- 1 ) 

100 Mev - 1 GeV 

1 - 10 GeV 

10- 100 GeV 

0017+8135 

S5 0014+813 

<0.33 

<3.0 x 10~ 3 


0035+5950 

1ES 0033+595 


<0.32 

<5.5 X 10^ 2 

0048+3157 

Mkn 348 

<0.49 



0120-2701 

lJy 01 18-272 

<0.70 

(1. 0+0.4) x 10~ 2 


0204-1701 

PKS 0202-17 

0.48+0.21 



0210-5101 

PKS 0208-512 

<0.64 

<1.3 x 10~ 3 


0214+5144 

GB6J 0214+5145 

<4.8 



0217+7349 

lJy 0212+735 




0217+0144 

PKS 0215+015 

0.91+0.31 

<7.7 x 10~ 3 


0221+3556 

lJy 0218+357 

<1.4 

<1.1 x 10~ 2 


0237+2848 

4C 28.07 

1. 5+0.5 

<1.1 x 10~ 2 


0238+1636 

PKS 0235+164 

<1.9 

<7.6 x 10^ 3 


0336+3218 

NRAO 140 

<0.95 

<1.5 x 10~ 2 

<4.5 x 10~ 4 

0423-0120 

PKS 0420-01 

1. 5+0.5 

<1.1 x 10^ 2 


0433+0521 

3C120 

<1.3 



0539-2839 

lJy 0537-286 

<0.59 



0738+1742 

PKS 0735+17 

<1.1 

(5.2+3. 1) x 10- 3 


0746+2549 

B2.2 0743+25 


<0.46 


0818+4222 

S4 08 14+425 

<1.2 

<1.4 x 10~ 2 


0824+5552 

OJ 535 

<0.39 



0841+7053 

4C 71.07 

<0.43 



0854+2006 

PKS 0851+202 

<0.81 

<1.1 x 10~ 2 


0915+2933 

B2 0912+29 

<0.46 

<7.8 x 10~ 3 


0923-2135 

PKS 0921-213 

<0.21 



0957+5522 

4C 55.17 

0.89+0.27 

(1. 7+0.5) x 10~ 3 

(l.l+l.l)x 10^ 4 

1058+5628 

1RXSJ 105837.5+562816 

0.51+0.19 

(6. 2+2. 3) x 10~ 3 


1058+0133 

4C 01.28 

<0.71 

<7.5 x 10~ 3 


1058-8003 

PKS 1057-79 

<0.65 

<7.4 x 10~ 3 

<9.7 x 10^ 5 

1131+3114 

B2 1128+31 

<0.24 



1136+7009 

S5 1133+704 

<0.77 

<8.0 x 10~ 3 


1153+4931 

4C 49.22 

<0.48 



1159+2914 

4C 29.45 

17.8+0.3 

(2. 0+0.5) x 10" 2 


1220+0203 

PKS 1217+02 

<0.50 



1221+2813 

ON 231 

<0.65 

<5.7 x 10~ 3 


1222+0413 

PKS 1219+04 

1. 2+0.4 

<7.3 x 10~ 3 


1246-2547 

PKS 1244-255 

<1.8 

<1.7 x 10^ 2 


1305-1033 

lJy 1302-102 

<0.55 



1310+3220 

lJy 1308+326 

<0.31 

(6. 5+2.2) x lO" 3 


1350+0940 

GB6B 1347+0955 

<0.43 



1419+0628 

3C 298.0 

<0.66 

<7.2 x 10~ 3 


1507+0415 

BZQJ 1507+0415 

<0.45 



1510-0543 

4C -05.64 

(45.5+7.6) x 10~ 2 

(21.3+2.9) x 10' 4 

<5.8 x 10~ 6 

1517-2422 

APLib 

1.22+0.41 

<9.5 x 10~ 3 

<4.8 x 10~ 4 

1555+1111 

PG 1553+113 

<1.1 

(14.6+3.3) x 10~ 3 

<7.8 x lO^ 5 

1625-2527 

OS -237.8 

<2.9 

<1.4 x 10~ 2 


1640+3946 

NRAO 5 12 

<9.8 

<0.67 

(5. 7+1. 7) x 10~ 2 

1719+4858 

ARP 102B 




1743+1935 

1ES 1741 + 196 

<0.31 


<2.5 x 10^ 2 

1840-7709 

PKS 1833-77 

<0.59 

<3.3 x 10~ 3 


1911-2006 

2E 1908.2-201 

<1.4 

<2.4 x 10~ 2 


1923-2104 

PMNJ 1923-2104 

1.97+0.67 

(8.78+5.4) x 10~ 3 


1924-2914 

OV -236 

<1.6 

(12.7+6.9) x 10~ 3 


2129-1538 

lJy 2126—158 


<6.1 x 10~ 3 

<4.2 x 10~ 4 

2147+0929 

lJy 2144+092 

1.47+0.45 

<7.6 x 10~ 3 


2151-3027 

PKS 2149-307 

<1.2 



2203+3145 

4C 31.63 

<0.32 



2207-5346 

PKS 2204-54 

<0.61 




Column 1 and 2: source name; Column 3-9: Fermi-LAT flux densities in 3 energy bands. 
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P. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 


Table 14. Fermi-LAT data at 1 bin/decade simultaneous with the Planck observations. Continued. 


J2000.0 name 

Source name 

Flux density (10 u 'phcm 2 s 

1 MeV- 1 ) 

100 Mev - 1 GeV 

1 - 10 GeV 

10- 100 GeV 

2209-4710 

NGC 7213 

<0.62 



2229-0832 

PKS 2227-08 

<0.81 

<9.7 x 10- 3 


2230-3942 

PKS 2227-399 

<0.31 


<9.1 X 10 5 

2303-1842 

PKS 2300-18 

<0.45 

<6.9 x 10 -3 


2327+0940 

PKS 2325+093 

<0.97 

<4.7 x 10 -3 


2333-2343 

PKS 2331-240 

<0.83 



2347+5142 

1ES 2344+514 


<9.9 x 10 -3 



Column 1 and 2: source name; Column 3-9: Fermi - LAT flux densities in 3 energy bands. 
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Table 15. Fermi-LK T data at 2 bins/decade integrated over 2 months around the Planck observations. 


P. Giommi et al.: Simultaneous Planck , Swift, and Fermi observations of blazars 
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Table 15. Fermi-LK T data at 2 bins/decade integrated over 2 months around the Planck observations. Continued. 


P. Giommi et al.: Simultaneous Planck , Swift, and Fermi observations of blazars 
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P. Giommi et al. : Simultaneous Planck, Swift, and Fermi observations of blazars 


Table 16. Fenni-LAT data at 1 bin/decade integrated over 2 months around the Planck observations. 


J2000.0 name 

Source name 

Flux density (10 lu phcm 2 s 

1 MeV~ 1 ) 

lOOMev - 1 GeV 

1 - 10 GeV 

10- 100 GeV 

0035+5950 

1ES 0033+595 

<0.41 

<3.8 x 10^ 3 

<7.0 x 10~ 5 

0048+3157 

Mkn 348 

<0.19 

<4.6 x 10^ 4 

<6.8 x 10^ 5 

0120-2701 

lJy 01 18-272 

(43.3+9.8) x 10~ 2 

(4. 7+ 1.2) x 10 -3 

(6.0+3.5) x 10- 5 

0214+5144 

GB6J 0214+5145 

<0.22 

<1.9 x 10^ 3 


0423-0120 

PKS 0420-01 

1.34+0.18 

(6. 5+1.5) x 10- 3 


0738+1742 

PKS 0735+17 

0.51+0.13 

(3.6+1.2) x 10 -3 

(7. 3+4.2) x 10- 5 

0854+2006 

PKS 0851+202 

0.50+0.14 

(3.2+1.19 x 10" 3 

(4. 1+2.9) x 10^ 5 

0915+2933 

B2 0912+29 

<0.18 

(96.5+5.4) x 10~ 5 


1305-1033 

lJy 1302-102 

<0.29 



1743+1935 

1ES 1741 + 196 

<0.39 

<2.9 x 10^ 3 


2347+5142 

1ES 2344+514 

<0.41 

<2.5 x 10^ 3 



Column 1 and 2: source name; Column 3-9: Fermi - LAT flux densities in 3 energy bands. 
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Table 17. Fermi-LA T data at 2 bins/decade integrated over 27 months. 


P. Giommi et al.: Simultaneous Planck , Swift, and Fermi observations of blazars 
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Table 17. Fermi-LA T data at 2 bins/decade integrated over 27 months. Continued. 


P. Giommi et al. : Simultaneous Planck, Swift, and Fermi observations of blazars 
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Table 18. Fenni-LAT data at 1 bin/decade integrated over 27 months. 


J2000.0 name 

Source name 

Flux density (10 lu phcm 2 s J 

1 MeV- 1 ) 

100 Mev - 1 GeV 

1 - 10 GeV 

10- 100 GeV 

0017+8135 

S5 0014+813 

<8.7 X 10 42 

<1.8 x 10~ 4 


0214+5144 

GB6J 0214+5145 

<6.7 x 10^ 2 

<4.2 x 10~ 4 

<6.4 X 10~ 6 

0217+7349 

ljy 0212+735 

(28.1+3.8) x 10~ 2 

(4.8+1.7) x 10' 4 

<3.7 x 10~ 6 

0550-3216 

PKS 0548-322 

<5.8 x 10^ 2 

<4.7 x 10~ 4 

(4.6+2. 8) x lO^ 6 

1147-3812 

PKS 1 144-379 

(23.1+3.0) x 10~ 2 

(15.7+2.2) x 10” 4 

<8.3 x 10~ 6 

1220+0203 

PKS 1217+02 

(12.3+3.6) x 10~ 2 

<4.9 x 10~ 4 

<5.8 x 10~ 6 

1305-1033 

ljy 1302-102 

<6.0 x 10^ 2 



1419+0628 

3C 298.0 

<0.10 

<4.0 x 10~ 4 


1423+5055 

BZQJ 1423+5055 

<4.2 x 10~ 2 



1456+5048 

1RXSJ 145603.4+504825 

<0.12 


<6.5 x 10~ 6 

1510-0543 

4C -05.64 

(45.5+7.6) x 10~ 2 

(21.3+2.9) x 10' 4 

<5.8 x 10~ 6 

1603+1554 

WE 1601 + 16W3 


<2.5 x 10“ 4 

<4.9 x 10~ 6 

1719+4858 

ARP 102B 

<4.8 x 10“ 2 

<1.4 x 10“ 4 


1743+1935 

1ES 1741 + 196 

<7.6 x 10~ 2 

(9. 5+1.9) x 10~ 4 

(5. 8+3.2) x 10’ 6 

2129-1538 

ljy 2126—158 

<0.11 

<3.5 x 10" 4 


2303-1842 

PKS 2300-18 


<2.7 x 10" 4 

<2.3 x 10~ 6 

2333-2343 

PKS 2331-240 

<9.5 x 10~ 2 

<4.3 x 10~ 4 

<3.7 x 10~ 6 


Column 1 and 2: source name; Column 3-9: Fermi - LAT flux densities in 3 energy bands. 
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